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Abstract. Symmetry group methods are applied to obtain all explicit group- invariant 



radial solutions to a class of semilinear Schrodinger equations in dimensions n ^= 1. Both 



focusing and defocusing cases of a power nonlinearity are considered, including the special 
^ ' case of the conformal power p = 4/n relevant for critical dynamics. The methods involve, 

firstly, reduction of the semilinear Schrodinger equations to group-invariant complex 2nd 
' order ODEs with respect to an optimal set of one-dimensional point symmetry groups, and 

secondly, use of inherited symmetries, hidden symmetries, and conditional symmetries to 
solve each ODE by quadratures. Through Noether's theorem, all conservation laws arising 
, from these point symmetry groups are listed. 

& . 

i. Introduction 

Semilinear Schrodinger equations with power nonlinearities for u(t,x), x G IR n , 

\u t = Am + k\u\ p u, Pt^O, k = const. 7^ (1.1) 

provide models of many interesting physical phenomena [T], such as propagation of laser 
beams in nonlinear media, slow modulation of plasma waves, motion of water waves at 
the free surface of an ideal fluid, dynamics of imperfect Bose condensates, and continuous- 
limits for mesoscopic molecular structure. In multi-dimensions n > 1, solutions exhibit 
very rich types of behaviour [21 [1] , particularly radial similarity solutions and more general 
group-invariant radial solutions which are important for investigating collapse or blow-up 
behaviour, dispersive behaviour, critical dynamics and asymptotic attractors, as well as for 
^ . testing numerical solution methods. 

To-date in the literature, only a few explicit n-dimensional radial solutions u(t, \x\) are 
apparently known (e.g. see Ref.[3]). Most of the work on explicit solutions to Schrodinger 
equations with power nonlinearities (II. ip has concentrated on systematically applying sym- 
metry methods [U E] to classify all possible types of group- invariant solutions [HI CO El E] , 
including radial and cylindrical as well as other less geometric types, only for the conformal 
power p = 4/3 in dimension n = 3 in addition to the lowest non-conformal even powers 
p = 2,4 in dimensions n = 2,3. These are the cases of greatest relevance for physical 
applications. 

The present paper, in contrast, will be devoted to deriving all group-invariant radial 
solutions u(t, \x\) of the Schrodinger equation ( 11. ip in all dimensions n^l and for all powers 
p 7^ 0, including the case of the general conformal power p = 4/n relevant for understanding 
critical dynamics. Both the focusing case k > and the defocusing case k < will be 
considered. 
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Group-invariant radial solutions for the Schrodinger equation (11. ip arise naturally from 
the corresponding radial Schrodinger equation 

vat = u rr + (n — l)r~ l u r + k\u\ p u (1-2) 

under reductions by one- dimensional groups of point transformations on r = \x\,t, u, u leav- 
ing the radial equation (jl.2p invariant. Each such symmetry reduction yields a complex 2nd 
order semilinear ODE 

U" = f(£,U,U,U',U') (1.3) 

formulated in terms of the invariants 

£ = E(t, r, u, u), U = Y(t, r, u, u), U = T(t,r,u,u) (1.4) 

determined by a given symmetry transformation group, provided that this system (jl.4p can 
be inverted (at least implicitly) to obtain both of the dependent variables u, u, and one of the 
independent variables r, t in terms of £, U, U, and the other independent variable (which will 
be the case whenever the orbits of the transformation group acting on the variables (r, t, u, u) 
are one- dimensional and have a projectable regular action on (r, £)). Then each solution of the 
ODE (11.31) for t/(£) will yield a group-invariant solution of the radial Schrodinger equation 
(TOD for u(t,r). 

If two groups of point symmetry transformations are related by conjugation with respect 
to some point transformation in the full symmetry group of the radial Schrodinger equation 
(11. 2p . then the action of this point transformation on solutions u(t,r) will map the group- 
invariant solutions determined by the two symmetry groups into each other. Consequently, 
for the purpose of finding all group-invariant radial solutions, it is sufficient to work with 
a maximal set of one- dimensional point symmetry groups that are conjugacy inequivalent. 
For each such group, once the resulting group-invariant solutions have been found, the full 
symmetry group of the radial Schrodinger equation can be applied on these solutions to 
obtain the group-invariant solutions determined by all other one-dimensional point symmetry 
groups in the same conjugacy class. 

For a given point symmetry reduction of the radial Schrodinger equation (ll.2p . the complex 
2nd order semilinear ODE ( II. 3p has an equivalent polar form given by a semilinear system 
of 2nd order real ODEs for the amplitude and phase of the variable 

U = Aexp(i$). (1.5) 

It is a non-trivial task to solve this semilinear system for A(£), $(£). One systematic approach 
is the method of reduction of order [U which relies essentially on the structure of the 
group of point symmetries admitted by the system and on the whether the system admits 
a Lagrangian structure or not, depending on the values of the nonlinearity power p and the 
dimension n. Firstly, equation ( ll.2p for u together with the complex-conjugate equation 
for u are the respective Euler-Lagrange equations SL/Su = and 5L/5u = of the radial 
Schrodinger Lagrangian 

L = (\xuu t - \m t u + u r u r - ^-ku 1+p / 2 u 1+p/2 )r n -\ (1.6) 

Thus, if this Lagrangian is invariant under the given point symmetry transformations used 
for the reduction, then the semilinear system for A(£) and $(£) wm inherit a Lagrangian 
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obtained by reduction of the radial Euler-Lagrange structure. Secondly, the phase rotation 
transformations 

u -»■ e 1 ^, u e-^u (1.7) 

comprise a !7(1) group of point symmetries of equation (ll.2p . This group is readily shown 
to commute with the full group of point symmetries admitted by equation (11.21) . Therefore, 
if the given point symmetry reduction starts from any other one-dimensional group of point 
symmetries, then the semilinear system for A(£) and will inherit a U(l) group of 

point symmetry transformations on £, A, $ obtained by reduction of the phase rotation 
transformations on u, u. Moreover, this system will also inherit all point symmetries that 
belong to the normalizer subgroup of the given point symmetry used for the reduction. An 
important remark is that the full group of point symmetries admitted by the system, or 
equivalently by the complex semilinear ODE for £/(£), can possibly contain "hidden" point 
symmetries in addition to those point symmetries inherited via reduction. 

In the case when a given point symmetry reduction of the radial Schrodinger equation 
(II. 2p inherits a Lagrangian structure, the semilinear system for A(£) and or equiva- 

lently the complex semilinear ODE for £/(£), can be solved explicitly (up to quadratures) 
if its Lagrangian is invariant under a two-dimensional group of point symmetry transfor- 
mations. In the alternative case when no Lagrangian structure exists, a four- dimensional 
point symmetry group with a solvable Lie algebra structure is needed to obtain the explicit 
solutions (up to quadratures) for A(£) and $(£)> or equivalently for t/(£)- I n all cases, each 
one- dimensional group of point symmetry transformations admitted by the semilinear sys- 
tem for A(£) and $(£)> or equivalently by the complex semilinear ODE for £/(£), can be used 
to find a single solution that is invariant under the admitted symmetry transformations on 
£, A, $, or equivalently on £, U, U. 

A different reduction method is applicable in the Lagrangian cases. If a given semilinear 
system for A(£) and inherits a [/(l)-invariant Lagrangian obtained by reduction of both 
the radial Schrodinger Lagrangian (11.61) and the U(l) group of phase rotation symmetries 
(II. 7p . then the 2nd order ODE for can be explicitly integrated so that the system 

reduces to a 2nd order ODE for A(£) alone, with $(£) given by an integral in terms of A(£) 
containing an arbitrary constant. This reduced ODE for A(£) has two useful features. First, 
nonlinear terms in the reduced ODE vanish if the nonlinearity power is p = — 1 and the 
arbitrary constant in integral for <&(£) is zero. In this case, the resulting linear ODE for 
can be solved explicitly in terms of special functions. Second, a Lagrangian structure can 
be derived for the reduced ODE starting from the [/(l)-invariant Lagrangian of the original 
semilinear system. Thus, in the nonlinear case this ODE can be solved (up to quadratures) if 
its Lagrangian is invariant under a one-dimensional group of point symmetry transformations 
on £, A. In all cases the reduced ODE also can be solved (up to quadratures) if it admits a 
two-dimensional group of non-variational point symmetry transformations. More generally, 
any one- dimensional group of point symmetries admitted by the reduced ODE can be used to 
find a single invariant solution for A(£). Each solution found for A(£) determines a solution 
for $(£) and hence yields a solution to the complex semilinear ODE for U(£). These solutions 
will turn out to differ from those solutions obtained via the previous reduction method if the 
group of point symmetries used to obtain A(£) is a "hidden" group which is not inherited 
under reduction from any point symmetries admitted by the semilinear system for A(£) and 
$(£), or equivalently by the complex semilinear ODE for £/(£). In particular, a "hidden" 
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symmetry group can arise from point transformations on £, A that leave invariant the integral 
for <H0- 

The rest of this paper is organized as follows. 

Section H] contains some preliminaries on symmetries and Noether's theorem. We first state 
the full point symmetry structure for both the n-dimensional Schrodinger equation (11.11) and 
the radial Schrodinger equation (11.21) . In addition we summarize the conservation laws that 
arise from the Lagrangian structure of these two equations via Noether's theorem. Next 
for the radial Schrodinger equation (jl.2p we present a maximal set of one-dimensional point 
symmetry groups that are conjugacy inequivalent. Modulo phase rotations, the symmetry 
groups in this set consist of time-translations, scalings, and a particular combination of 
inversions (conformal transformations) and time-translations. These three symmetry groups 
will be referred to as the optimal subgroups for symmetry considerations. We then write 
down the complex 2nd order semilinear ODEs (11.31) given by reduction under each point 
symmetry subgroup in the optimal set, and we summarize the Lagrangian structure admitted 
by each of the ODEs, depending on p and n. 

Section [3] explains the reduction of order method in detail for solving complex 2nd order 
semilinear ODEs (11.31) by use of point symmetries. In particular, we show how to streamline 
the standard reduction steps in an efficient way by combining the Lagrangian and non- 
Lagrangian cases through the use of canonical coordinates determined by any admitted 
one-dimensional point symmetry group. We also explain details of the alternative reduction 
method in the £7(l)-invariant Lagrangian cases, which we carry out by again using canonical 
coordinates to streamline the steps. 

Sections HJ \5\ M apply these reduction methods to derive solutions U(£) for each ODE 
(II. 3p arising from the three optimal subgroups of point symmetries for the radial Schrodinger 
equation (jl.2p . In particular, we are able to obtain explicit solutions in terms of elementary 
functions. 

The main results of the paper are presented in section [7J We first list all of the group- 
invariant radial solutions u(t, r) determined by applying the full group of point symmetries 
of the radial Schrodinger equation ( II .2p to each of the group-invariant solutions U(£) derived 
from the three optimal subgroups of point symmetries. We next discuss the analytical fea- 
tures of these solutions u(t,r). Some of the solutions and underlying symmetries are found 
to involve non-integer values of n. We show that for such cases the radial Schrodinger equa- 
tion can be interpreted alternatively as modeling two-dimensional radial wave propagation 
with a point-source disturbance acting at the origin, where the strength of the disturbance 
is proportional to the value of n — 2. 

Finally, concluding remarks are made in Section [HJ 

2. Symmetries and Conservation Laws 

For the Schrodinger equation ( II .ip in M™, a point symmetry is a one-dimensional Lie 
group of transformations acting on the variables (t, x, u, u) such that the prolongation of its 
infinitesimal generator 

X = r(t, x, u, u)d/dt + ((t, x, u, u) ■ d/dx + r](t, x, u, u)d/du + fj(t, x, u, u)d/du (2.1) 

satisfies prX(ru t — Am — k\u\ p u) = for all formal solutions u(t,x) of equation ( II. ip . This 
is the condition for equation (II. ip to be infinitesimally invariant under the transformation 
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group generated by X. Each such generator acting on solutions u(t,x) has an equivalent 

characteristic form 



X = P(t, x, u, u, u t , u t , Vw, Vu)d/du + P(t, x, u, u, u t , u t , Vtt, Vu)d/du (2.2) 

with 

p = n - ru t - c • Vw (2.3) 

satisfying 

iD t P - D x ■ D X P - k(l + p/2)u p/2 u p/2 P - k(p/2)u p/2 ~ 1 u 1+p/2 P = (2.4) 

as given by applying prX to equation (11 .11) and then eliminating ut, Ut, and x-derivatives of 
u t , u t through the equation (II. ip on u(t,x) and the complex-conjugate equation on u(t,x). 
Here D t and D x denote total derivatives with respect to t and x. Thus, a solution u = f(t, x) 
of the Schrodinger equation (11.11) is group-invariant under a one-dimensional point symmetry 
group with a generator (12.11) if (and only if) it satisfies the additional equation 

r](t, x, u, u) — r(t, x, u, u)u t — C(^j x , u , u) ■ Vw = 0. (2.5) 

It is straightforward to solve equation (12. 4p to determine all point symmetry generators 
(12. ip . In particular, because r, (, r] do not depend on any derivatives of u and u, the 
equation (12. 4p splits with respect to x-derivatives of u and u, yielding an overdetermined 
linear system of PDEs on the functions r(t,x,u,u), ((t,x,u,u), r](t,x,u,u). This system 
yields the following well-known result [6l [10] . 

Theorem 1. The point symmetries of the Schrodinger equation (11. ip are generated by 

phase rotation Xi = md/du — md/du, (2.6a) 

time translation X 2 = d/dt, (2.6b) 

space translations ^~3(i) = ■ d/dx, I — 1, ... ,n, (2.6c) 

Galilean boosts X 4 (^ = 2te^ ■ d/dx — i(e^) • x){ud/du — ud/du), I = 1, . . . ,n, 

(2.6d) 

space rotations X 5 (^/) = (e^) ■ x)e^') ■ d/dx — (e^ ■ x)e^ ■ d/dx, 

1 = l,...,n-l;l' = l + l,...n, (2.6e) 
scaling X 6 = 2td/dt + x ■ d/dx — (2/p)ud/du — (2/p)ud/du, (2.6f) 
inversion X7 = t 2 d/dt + tx ■ d/dx — (2t/p + i\x\ 2 / A)ud / du 

— (2t / p — i\x\ 2 / A)ud / du only for p = A/n, (2.6g) 



where {em, . . . , et n )} is any orthonormal basis forW" 1 . The corresponding transformation 
groups acting on solutions u = f{t,x) of the Schrddinger equation (II. ip are given by 

u = exp (i0)/(t, x), (2.7a) 

u = f(t-e,x), (2.7b) 

u = f(t,x-ee(i)), 1 = 1, ...,n, (2.7c) 

u = exp (— i(ee(j) • x/2 - e 2 t/4)) f(t,x - ete^), l = l,...,n, (2.7d) 

u = f(t, x + (cos(0) - l)(e(j) e (z) + e (P) e (F) ) • x + 2 sin(</>)(e w A e^) • x), 

I = l,...,n-l;l' = l + l,...n, (2.7e) 

u = \- 2 / p f(\~ 2 t,\- l x), (2.7f) 

u = (1 + ety 2/p exp (-ie|a;| 2 /(4 + Aet)) /(t/(l + et), x/(l + et)) only for p = A/n, 

(2-7g) 

with group parameters — oo < e < oo ; < A < oo, < < 2tt. In the case p ^ A/n, the 
transformations (I2.7al) - (12.7fp comprise a semi-direct product of a scaling group acting on a 
central extension of the Galilean group, which has the Lie algebra structure 

[Xi,X 2 ] = [Xi,X 3 (i)] = [Xi,X 4 (j)] = pCi,X 5 (j )J /)] = 0, (2.8) 

[X2,X 3( z)] = 0, [X 2 ,X 4 (j)] = 2X 3 (i), [X 2 ,X 5 ( ii ^)] = 0, (2.9) 

[X3(fc), X 4 (;)] = SkiXi, [X 3 (fc), X 5 (j ) p)] = 5kiX 3 (i/) — 5fci'X 3 ( Z ), [X 4 (fc), X 5 (/ ) //)] = 5ki^i(v) — Sm'X^i), 

(2.10) 

[X 6 ,X!] = 0, [X 6 ,X 2 ] = — 2X 2 , [X 6 ,X 3 ( fc )] = — X 3 ( fc ), [X 6 ,X 4 ( fc )] = X 4 ( fc ), [X 6 ,X 5 (^/)] = 0. 

(2.11) 

In the case p = A/n, the inversion transformation ( |2.7g ) intertwines non-trivially with the 
previous group, as given by the commutator structure 

[X 7 , X x ] = [X 7 , X 4(0 ] = [X 7 , X 5(M0 ] = 0, (2.12) 

[X 7 ,X 2 ] = — X 6 , [X 7) X 3 (j)] = -|X 4(/ ), [X 7 ,X 6 ] = — 2X 7 . (2.13) 

Here © and A respectively denote the symmetric and antisymmetric parts of an outer 
product. The special power p = A/n for which the inversion group exists is commonly called 
the conformal power. 

For any domain Q C R™, the Schrddinger equation ( ll.ip has the variational formulation 

I - < 2 - 14 > 

given by the Lagrangian functional 

£= / / L{u,u,u t ,u t ,Vu,Vu) d n xdt, L = ~\uu t — \mu t + | Vw| 2 |^A;|u| 2+p . (2.15) 

A variational point symmetry of this functional (12. 15ft is an infinitesimal point transformation 
(12. ip on (t, x, u, u) under which £ is invariant up to spatial boundary terms at dQ and 
temporal boundary terms at t = to and t = t\. This invariance condition holds if and only 
if the Lagrangian satisfies 

prXL = D t A + D x ■ B (2.16) 
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for some functions A and B of t, x, u, u, and derivatives of u, u with respect to t and x, 
where X is the characteristic form fl2.2p - fl2.3l) of the generator (12.1 p . An equivalent condition 
on the Lagrangian is that prXL is annihilated by the variational derivatives with respect to 
u and u. 

Since invariance of C implies that its extrema (12. 14j) are preserved, every variational point 
symmetry of the Lagrangian functional (I2.15P for the Schrodinger equation (II -ip is thereby 
a point symmetry of the Schrodinger equation itself such that 

<5(prXL) <5(prXL) 
V %_ ; =0, J - = 0. 2.17 

ou ou 

This provides a straightforward way to determine all of the variational point symmetries 
starting from Theorem [H which yields the following result. 

Theorem 2. The variational point symmetries of the Schrodinger equation (II. ip are gen- 
erated by phase rotation f )2.6ap . time translation (12.6bp . space translations (I2.6c|) . Galilean 
boosts fl2.6dp . space rotations (I2.6e[) . scaling fl2.6fp only for the conformal power p = 4/n, 
and inversion ( |2.6gP for the conformal power p — 4/n. 

Variational symmetries give rise to conservation laws for the Schrodinger equation fll.ip 
by means of Noether's theorem as follows. The invariance condition (I2.16P combined with 
the variational identity 

5L „ 5L - „ / dL „ dL ^\ „ / dL „ dL - 



V 1 ^-^ fi u ^ ' ' ~*~ ^ >t \Q Ut ^' ' ^ du t ^y ^ x \dVu^ dVu / (2.18) 
yields the multiplier equation 

Re((m t - Ah - A;|u| p m)P) = AT + A • X (2.19) 

with 

r = X-^P-|£p, X = B -^P-JLP. (2.20) 

OU t OU t o\u OVU 

On all formal solutions u(t,x) of the Schrodinger equation fll.ip . the multiplier equation 
f)2.19p then produces a conservation law 

D t T + D x -X = 0, (2.21) 

where the conserved density T and flux X are given by the Noether relation (12.201) . In 
particular, A and B can be shown to have the specific form 

A = -tL, B = -(L (2.22) 

as derived from the invariance condition f)2.16p . Therefore, the conserved density and flux 
are given explicitly by the simple formulas 

T = —tL + Re(mP), X = -(L - 2Re(PVw) (2.23) 

in terms of L and P = r] — ru t — ( ■ Vu. This result together with Theorem [2] yields all 
of the conservations laws generated by the variational point symmetries of the Schrodinger 
equation (11.11) . as shown in Table [D using the notation V(z) = ■ V where {e(i), . . . , e( n )} is 
any orthonormal basis for R™. 
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Table 1. Conservation laws of the semilinear Schrodinger equation (II. ip in W 



2.1. Radial reduction. A space rotation generator (I2.6ej) acts as an infinitesimal rotation 
in the 2-plane specified by a pair of basis vectors em and em) in M n with I ^ V . The 
corresponding group of transformations 

(t, x, u) (t, x + (cos(0) - l)(em em + em) em)) • x + 2 sin(0)(e w A em)) ■ x, u) (2.24) 

is an SO(2) Lie group. Composition of all such transformations acting in the n(n — l)/2 
distinct 2-planes determined by an orthonormal basis {e(i) 3 • • • , e( n )} of M. n produces an 
SO(n) Lie group of rotations, whose invariants are functions of t, \x\ = r, u, u. 

Reduction of the Schrodinger equation (11.11) in M" under this point symmetry group of 
rotations gives the radial Schrodinger equation (11.21) where u(t,r) is a group-invariant so- 
lution. The radial Lagrangian (11.61) arises naturally from this reduction of the Lagrangian 
functional (I2.15P in W 1 , as given by 

rh roc 

£ ra d. = / / L(r,u,u,u t ,ut,u r ,u r ) drdt (2.25) 
J t Jo 
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in the radial domain < r < oo. 

From general results on symmetry reductions [I], a point symmetry generator (12. ip of the 
Schrodinger equation (11.11) in IR n will admit a radial reduction if and only if it belongs to the 
normalizer algebra of the so(n) Lie subalgebra of space rotations (I2.6ep in the Lie algebra of 
all point symmetry generators (12. 6p . Likewise, through Noether's theorem, a conservation 
law (I2.2ip of the Schrodinger equation ( II. ip will admit a radial reduction if and only if its 
corresponding variational point symmetry generator belongs to this same normalizer algebra. 

Because the set of all formal radial solutions u(t, r) is contained as a subset in the set of 
all formal solutions u(t,x) to the Schrodinger equation (II. ip in R n , the radial Schrodinger 
equation ( II .2p could possibly admit additional point symmetries and conservation laws other 
than those inherited through radial reduction. Consequently, the only way to find all radial 
point symmetries, as well as all radial conservation laws corresponding to variational point 
symmetries, is by directly solving the determining equations for their generators. 

For the radial Schrodinger equation (ll.2p . the generator of a point symmetry acting on 
the variables (t, r, u, u) is given by 

X = r(t, r, u, u)d/dt + p(t, r, w, u)d/dr + r/(t, r, w, u)d/du + fj(t, r, w, u)d/du (2.26) 

such that prX(m t — u rr — (n — l)r~ l u T — k\u\ p u) = holds for all formal solutions u(t, r) of 
equation (ll.2p . The characteristic form for each such generator 

X = P(t, r, u, it, u t ,u t ,u r , u r )d/ du + P(t, r, u, u, u t ,u t ,u r , u r )d/ du (2.27) 

with 

P = 7] - TUt - pUr (2.28) 

satisfies 

iD t P - D r 2 P - (n - ly^DrP - k(l + p/2)u p/2 u p/2 P - k(p/2)u p/2 - 1 u 1+p/2 P = 0, (2.29) 

where u t , Ut, and r-derivatives of u t , u t are eliminated through the equation (II. 2p on u(t, r) 
and the complex-conjugate equation on u(t,r). Here D t and D r denote total derivatives 
with respect to t and r. This determining equation (12.291) splits with respect to r-derivatives 
of u and u, giving an overdetermined linear system of PDEs on the functions r(t,r,u,u), 
p(t,r,u,u), r](t,r,u,u). After a straightforward integrability analysis, the system reduces to 
the PDEs 



T r = r u = T Q = p u = pu = r nu = Vu = 0, (2.30a) 

rjuu = Van = (pn - A)t u = T tu = 0, (2.30b) 

P = \rr u Vut = Vat, Vur ~ Vur = ~iPt, Vu + Vu = ~~Pr, (2.30c) 

ur] + ufj = uu(r] u + fj u ) , (2.30d) 

which are easily solved to obtain the following result. 
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Theorem 3. The point symmetries of the radial Schrodinger equation (jl.2p are generated 
by 

phase rotation X p h as . = md/du — md/du, (2.31a) 

time translation X trans . = d/dt, (2.31b) 

scaling X sca i. = 2td/dt + rd/dr — (2/p)ud/du — (2/p)ud/du, (2.31c) 
inversion Xi nvcr . = t 2 d/dt + trd/dr — (2t/p + ir 2 / A)ud / du 

— (2t/p — ir 2 /A)ud/du only for p = A/n, (2.31d) 

with the Lie algebra structure 

[Xphas. > Xtrans.] [Xphas. j X sca i ] [Xphas. ; Xj nv er.] 0, (2.32) 
[Xtrans. ) X sca l.] 2X^ ranSi , [X^ ranSi , Xj nv er.] X sca l. ; [X sca h , Xj nvcr J 2Xj nver . (2.33) 

The corresponding transformation groups acting on solutions u = f(t, r) are given by 

w = exp(i0)/(t,r) J (2.34a) 

u = f(t-e,r), (2.34b) 

u = \- 2lp f(\-H,\- l r), (2.34c) 

u = (l + et)- 2/p exp(-ier 2 /(4 + 4et))/(t/(l + et),r/(l + et)) only for p = A/n, (2.34d) 
with group parameters — oo < e < oo ; < A < oo ; < <p < 27r. 

Thus, the only point symmetries admitted by the radial Schrodinger equation ( 11. 2p are 
the ones it inherits from reduction of the Schrodinger equation ( II. 2p in M n under the SO(n) 
group of rotations. 

The Lagrangian functional (I2.25P for the radial Schrodinger equation ( II. 2p is invariant up 
to spatial boundary terms (at r = and r — > oo) and temporal boundary terms (at t = t 
and t = t\) under an infinitesimal point transformation ( I2.26P on (t, r, u, u) if and only if 

prXL = D t A + D r B (2.35) 

holds for some functions A and B of t, r, u, u, and derivatives of u, u with respect to t 
and r. This is the condition for X to be a variational point symmetry of the Lagrangian 
functional ( I2.25p . Since the transformation group generated by a variational point symmetry 
necessarily preserves the extrema SC^d./Su = of this functional £ ra d., every such symmetry 
generator X is thereby a point symmetry of the radial Schrodinger equation itself for which 
the variational condition (I2.35P holds, or equivalently for which the equation 

du du 

is satisfied. Condition ( I2.36P can be shown to reduce to the symmetry determining equations 
(I2.30p plus the additional equation 

r t + (n - l)r~V - (1 + 4/p)p r = 0. (2.37) 

It is straightforward to determine which of the point symmetries from Theorem [3] satisfy this 
equation (12.371) . which yields the following result. 
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Theorem 4. The variational point symmetries of the radial Schrddinger equation (jl.2p are 
generated by phase rotation ( 12. 3 lap . time translation (12.31b[) . scaling (I2.31cl) only for the 
conformal power p = 4/n, and inversion (12.31d|) /or the conformal power p = 4/n. 



From Noether's theorem, these variational symmetries yield the multiplier equation 



l)r 1 u r 



k\u\ p u)P) = D t T + D r X, 



(2.38) 



Ke((iu t — u„ — (n 
which produces a conservation law 

D t T + D r X = (2.39) 

holding on all formal solutions u(t, r) of the radial Schrddinger equation (ll.2p . Here P is the 
function (I2.28P given by the characteristic form of a symmetry generator (I2.27p . while the 
conserved density T and flux X are given by the simple formulas 

T = -rL + Re(wP), X = -pL - 2Re(Pw r ), (2.40) 

which can be derived from the variational condition ( I2.35P in terms of the radial Lagrangian 
(11. 6p . The resulting conservations laws generated by all of the variational point symmetries 
in Theorem S] are shown in Table [2J 



T 



-r n 1 u 2 



X 



L 2 norm 



lr 



- n_1 (kl a -^L%l^ 



(uU r — MM r ) 



-r n 1 (u t u r + u t u r ) 



energy 



r n (\Ur\ 2 + d2 k \ U \ P ) 



2tr n -\\u r \ 2 - £- 2 k\u\ 2+ *) 
+ \\r n {uu r — uu r ) 



hr n (uu t 



2„n-l 



p+2' 
UU t ) 



-r 



{uU r + MM, 



dilation energy 
p = 4/n 



2tr n {u t u r + u t u r ) 



— Hr n 1 {uu r + uu r ) 

— t 2 r n ~ l {u t u r + u t u r ) 
-tr n (\u r \ 2 + £- 2 k\u\ 2+ v) 
+ \itr n {uu t — UUt) 

_l_ I ir ra+1 (^ Ur _ U u r ) 



2 ""-W - £- 2 k\u\ 2+p ) 



t 2 r 



+ \itr n {uu r — uu r 
+ \r n+1 \u\ 2 



conformal energy 
p = 4/n 



Table 2. Conservation laws of the semilinear radial Schrddinger equation ( 11.2 



2.2. Group invariance of radial solutions. A solution u = f(t,r) of the radial 
Schrddinger equation (II. 2p is group-invariant under a one-dimensional point symmetry group 
if (and only if) it satisfies X.u\ u _^ t . — where X is the symmetry generator in characteristic 

form (I2.27I) - (I2.28I) . To find all group-invariant solutions, it is sufficient to consider a maximal 
set of one-dimensional point symmetry groups that are conjugacy inequivalent in the full Lie 
group of point symmetries admitted by the radial Schrddinger equation (11. 2p . This is most 
easily carried out at the Lie algebra level by exhibiting a maximal set of one- dimensional 
subalgebras that are conjugacy inequivalent in the Lie algebra of point symmetry generators. 
Such a set, called optimal symmetry generators, can be straightforwardly determined by the 
methods in Refs. [H [TT] starting from the point symmetry algebra given in Theorem [3J 
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Lemma 1. The Lie algebra of point symmetry generators (12.311) for the radial Schrddinger 
equation (11. 2p is isomorphic to u(l) © A(2) when p ^ 4/n and u(l) ©sl(2,R) when p = 
4/n. (Here A(2) denotes a 2- dimensional, non-abelian Lie algebra, which is unique up to 
isomorphism.) An optimal set of symmetry generators X consists of 

X phas ., (2.41a) 

Xtrans. + ^Xphas.j (2.41b) 
Xscai. + /iX phas ., (2.41c) 

with parameters — oo < v < oo, — oo < /x < oo ; and afeo 

Xtrans. Xi nvcr + KX phas . for p = A/n (2.41d) 

with parameter — oo < k < oo. 

These generators (12.411) define group- invariant solutions determined by the respective equa- 
tions 

u = 0, (2.42a) 
iz/ M - M t = 0, (2.42b) 
( 2 -ifx)u + 2tu t + ru r = 0, (2.42c) 
(2* + i(|r 2 - + (£ 2 + + iru r = 0, p = 4/n, (2.42d) 

which come from prXn = 0, where X denotes a generator (12.411) in characteristic form 
(I2.27p — (I2.28p . The corresponding point transformation groups are given by 

t — > t, r — > r, u — > exp(i0)n, (2.43a) 
t — > t + e,r — > r, u — > exp(iz/e)n, (2.43b) 

t -»■ A 2 t, r -> Ar, u -> exp(i/z In A) A~ 2/p n, (2.43c) 
t — > (sin + t cos 0) / (cos — t sin 0) , r — )■ r / (cos — t sin 0) , 

it — > exp (ik0 — i|r 2 sin0/ (cos0 — t sin 0)) (cos — i sin 0) 2 / p n, p = 4/n, (2.43d) 

with group parameters — oo <e<oo, 0<A<oo, < < 2n. 

Invariance under the phase rotation symmetry generator (j2.41aj) yields only a trivial group- 
invariant solution ( I2.42ap . Hence, for finding non-trivial group-invariant solutions, only the 
remaining symmetry generators fl2.41c[) — (12.41d[) in the optimal set need to be considered. 



2.3. Radial scaling reductions. Consider the symmetry group of scaling transformations 
(I2.43cj) with the generator (|2.41c[) . Its invariants (ll.4p satisfying X sca i.£ + /iX pn as.£ = 0, 
Xgcai.C/ + fjX phaa .U = and X sca i.C7 + /iX phas .f7 = are given by 

£ = t/r 2 , U = r 2/p exp(-i/ilnr)n, U = r 2/p exp(i/ilnr)w, (2.44) 

depending on a parameter — oo < u < oo. Hence 



u(t,r) = r 2//p exp(i/ilnr)?7(£) 
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(2.45) 



yields the corresponding form for group-invariant solutions. The radial Schrodinger equation 
(11. 2p thereby reduces to the complex semilinear 2nd order ODE 

4£ 2 U" + ((8 - 2n + 8/p)£ - i(l + A^))U' 

+ ((4 - 2n)/p + A/p 2 - /i 2 + in(n - 2 - A/p))U + k\U\ p U = 0. (2.46) 

For fi = 0, this reduction describes radial similarity solutions. 

From Theorem HJ the symmetry group (I2.43cl) leaves invariant the Lagrangian functional 
(I2.25P for the radial Schrodinger equation iff p = 4/n. Hence, in the case of the conformal 
power, ODE (12.461) has a variational formulation given by reduction of the radial Schrodinger 
Lagrangian (II. 6p . 

L = -4e\U'\ 2 + (^-^)\U\ 2 + ^ r2 k\U\ 2+ P + i^ + 2^)(UU'-UU') for p = A/n. (2.47) 

2.4. Radial time-translation reductions. Consider the symmetry group of time- 
translation transformations (12.43bl) with the generator (12.41bj) . Its invariants (11.41) satisfying 
Xtmnsi + ^X phas .£ = 0, X tvans U + uX phas U = and X trans U + uX phas U = are given by 

£ = r , U = exp (-wt)u, (2.48) 

depending on a parameter — oo < v < oo. Hence 

u(t,r) = exp(wt)U(£) (2.49) 

yields the corresponding form for group-invariant solutions. The radial Schrodinger equation 
(II .2p thereby reduces to the complex semilinear 2nd order ODE 

U" + (n- + uU + k\U\ p U = 0. (2.50) 

For v = 0, this reduction describes radial static solutions, which are invariant under the 
symmetry subgroup of time-translation transformations 

t — y t + e,r — y r, u — y u. (2-51) 

Since, by Theorem HI the Lagrangian functional (12.251) for the radial Schrodinger equation 
is invariant under the symmetry group (I2.43b|) . ODE (I2.50P has a variational formulation 
given by reduction of the radial Schrodinger Lagrangian (11.61) . 

L = e~ l {-\U'\ 2 + v\U\ 2 + ^- 2 k\U\ 2+p ). (2.52) 

2.5. Radial conformal reductions. Consider the symmetry group of combined inversion 
and time-translation transformations ( I2.43dp with the generator (I2.41d.fl for p = A/n. Its 
invariants (JIIJ) satisfying X trans £ + X im/CT £ + nX phas ^ = 0, X trans U + X irwer U + nX phas .U = 
and X tvans U + X inver .[/ + nX phas U = are given by 

£ = (1 + t 2 )/r 2 , U = r 2/p exp (ik arctan(l/t) + ir 2 t/(4(l + t 2 )))u, 

U = r 2/p exp ( - i«arctan(l/t) - ir 2 t/(4(l + t 2 )))u, 

depending on a parameter — oo < k < oo. Hence 

u(t, r) = r~ 2/p exp ( - in arctan(l/t) - ir 2 t/(4(l + t 2 ))) U(£) (2.54) 

yields the corresponding form for group-invariant solutions. The radial Schrodinger equation 
(II. 2p thereby reduces to the complex semilinear 2nd order ODE 

Ai 2 XJ" + 8£U' + (kC 1 - r 2 /4 + n(l - n/A))U + k\U\ 4/n U = 0. (2.55) 
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From TheoremHJ since the Lagrangian functional f!2.25j) for the radial Schrodinger equation 
is invariant under the symmetry group (12.43dj) . ODE (I2.55P has a variational formulation 



given by reduction of the radial Schrodinger Lagrangian (11.61) with p = 4/n, 

L = -±e\U'\ 2 + {kC 1 - IC 2 + n(l - \n))\U\ 2 + ^k\U\ 2+ " n ). (2.56) 

2.6. Group invariance of optimal radial reductions. The ODEs (I2.46p . (I2.50p . (I2.55P 
arising by the reduction of the radial Schrodinger equation under its optimal point symmetry 
subgroups (I2.41cl) . (12.41bl) . f!2.41dp have the following symmetry structure. 

Proposition 1. The only point symmetries admitted by the scaling-group ODE (I2.46P and 
the conformal- group ODE (I2.55P are generated by phase rotations 

Y p has. = lUd/dU - lUd/dU. (2.57) 

The point symmetries admitted by the translation-invariant ODE (12.501) are generated by 

scaling Y scal . = £d/d£ - (2/p)Ud/dU - (2/p)Ud/dU for v = 0, (2.58) 

dilation Ydii. = £ 2 ' r \£/{n - 2)d/d£ - Ud/dU - Ud/dU) (2.59) 

for v = 0, p = 2(3 - n)/(n - 2) ^ 0, 

in addition to phase rotations (12.571) . 

In the case v = 0, the translation-group ODE (12.501) determines all static solutions of the 
radial Schrodinger equation (II. 2p . 

U" + (n - l^U' + k\U\ p U = with U = u,£ = r. (2.60) 

The scaling symmetry (12.581) of this ODE is inherited from the scaling invariance of the 
radial Schrodinger equation, due to the commutator structure [X trans . , X scal J = 2X trans .. In 
contrast, the dilation symmetry (12.591) is not inherited from any invariance of the radial 
Schrodinger equation and thus describes a hidden symmetry arising for static solutions. 

3. Quadrature of complex 2nd order semilinear ODEs 

Details of the two reduction of order methods outlined in section [1] for solving complex 
2nd order semilinear ODEs of the general form (12.461) . (I2.50p . (12.551) will now be presented. 

3.1. Reduction by point symmetries. Consider a complex 2nd order semilinear ODE 

<x(Z)U" + + l(0U + k\U\?U = (3.1) 

with a real independent variable £ and a complex dependent variable U . Here a, (3, 7 are 
allowed to be complex functions of £, while k and p are assumed to be non-zero real constants. 
Every such ODE (13. ip is invariant under a U(l) group of phase rotation symmetries given 
by the generator (12.571) . Now suppose an ODE (13.11) admits another one-dimensional group 
of point symmetries, with a generator of the form 

Y = C(09/dC + n(£)Ud/dU + Q(0Ud/dU, (3.2) 

whereby 

[Y phas .,Y] = 0. (3.3) 
If the ODE has a variational formulation for which both Y p h as and Y are variational sym- 
metries, then these two symmetries can be used to reduce the ODE to quadratures by means 
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of first integrals. However, the standard reduction steps |4J would require finding the La- 
grangian functional and checking its invariance under the two symmetries, which can be 
cumbersome to carry out. A simpler, more direct way to accomplish the same reduction is 
to utilize canonical coordinates and integrating factors associated to the symmetries Y phas . 
and Y as follows. 

Make a change of variables from (£,£/, U) to (z, V, V) given by the canonical coordinates 
of the generator (13. 2p . 

Yz = 1, YV — 0, YV = 0, (3.4) 

where 

z = J (1/CR, V/U = exp ( - J (Sl/Odt) (3.5) 

are functions only of £. In terms of these variables, the generator (13. 2p takes the form of a 
z-translation 

Y = d/dz, (3.6) 

while the phase rotation generator (I2.57P is given by 

Y phas . = iVd/dV - iVd/dV. (3.7) 

Consequently, the ODE (13. ip is transformed into 

aV" + bV + cV + k\V\ p V = (3.8) 

for V(z), where the coefficients a,b,c are constants. The following result is straightforward 
to prove by considering the characteristic form of the symmetry generators (I3.6p - (l3.7p . 

Y = -V'd/dV - V'd/dV, Y phas . = iVd/dV - iVd/dV. (3.9) 

Lemma 2. (i) An ODE (13.81) has a non-trivial invariant solution with respect to z- 
translations iff c — c and, when p is even, c/k < 0. Then V — yields 

\V\ = (-c/k) 1 ^, (3.10) 

which determines V up to an arbitrary constant phase, (ii) An ODE (13. 8p has integrating 
factors V' and —lV iff a = a, b = —b, c = c. Then the first integrals are respectively given 
by 

a\V'\ 2 + c\V\ 2 + 2k J \V\ p+1 d\V\ = Ci, (3.11) 

a(VV - V'V) + b\V\ 2 = i2C 2 , (3.12) 

from which V(z) can be determined by quadratures, where C\ and Ci are arbitrary real 
constants. 

The solutions arising from Lemma [2] can be easily obtained in a more explicit form in 
terms of polar variables 

V = Aexp(i^). (3.13) 
In particular, the first integrals (13 . 1 1 [) — (13. 12[) yield 

aA' 2 = -aA 2 & 2 - cA 2 -2k J A p+1 dA + C u (3.14) 

a$' = b + C 2 /A 2 , b = i2b. (3.15) 
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Combining these two differential equations, we get the quadratures 

dA 



±z + C 3 , (3.16a) 
dA 



v^P) 

/dA 
- + C 4 , (3.16b) 
a 2 \/hTa) 



assuming A' ^ 0, where 

H(A) = (ad - 2bC 2 )/a 2 - ((b/af + c/a)A 2 - (C 2 /a) 2 A~ 2 - (2k/ a) I A p+1 dA. (3.17) 



The quadratures (I3.16P determine the general solution V(z) of ODE (13. 8p in the case a = a, 
b = —b, c = c. Moreover, these conditions on a, b, c are necessary and sufficient for the ODE 
(13. 8p to have an Euler-Lagrange structure 

5L/5V = aV" + bV + cV + k\V\ p V = 0, a = a,b = -b 7 c = c (3.18) 

given by the Lagrangian 

L = -aV'V + \b(VV' - V'V) + \c\V\ 2 + 2k J \V\ p+1 d\V\, (3.19) 

for which both the ^-translation generator (13. 6p and the phase rotation generator (13 .7p are 
variational symmetries. Noether's theorem thereby asserts 

prYL = -f z = -V> 5 ± - V' b ± + ±((aV> - \bV)V + (aV + \bV)V% (3.20) 

XT XT A 

prY phas .L = = iV— - iV-^ + -^{(-aV' + \bV)\V + (aV + \bV)\V), (3.21) 

giving a derivation of the first integrals (I3.1ip -( l3.12p . 

In the general case without conditions on a, b, c, the ^-translation invariant solution V = 
const, of ODE (13. 8p is given by 

A = (-c/k) 1/p = const., $ = const.. (3.22) 

Finally, changing variables (z, V, V) back into (£, U, U) in the original complex semilinear 
ODE (13. ip yields (i) the group-invariant solution with respect to the symmetry generator 
( 13.2 j) ; (ii) the general solution under necessary and sufficient conditions for the generator 
(13. 2 p to be a variational symmetry. 

3.2. Reduction by hidden conditional symmetries. Consider a real 1st order La- 
grangian of the form 

L = w(0{-HOUV' + 1 m(U'U-U'U)+m\U\ 2 + 2k I \U\ p+1 d\U\), (3.23) 

where zu, a, (3, 7 are real functions of £, and k 7^ 0, p 7^ are real constants. The Euler- 
Lagrange equation SL/SU = yields a complex 2nd order semilinear ODE of the form (13. ip 
whose coefficients a, (3, 7 are given by 

a = a, = a' + aw'/w - i2/3, 7 = 7 - 10' + (3w'/w). (3.24) 

16 



Both the Lagrangian (I3.23j) and the Euler-Lagrange ODE (13. ip . (I3.24p are invariant under 
the generator (I2.57P of phase rotation symmetries on U, U . Hence Y p h as . is a variational 
symmetry to which Noether's theorem can be applied to obtain a reduction of this ODE. 
The simplest way to carry out the reduction is by use of polar variables 

C/ = Aexp(i$). (3.25) 

In particular, the ODE (13. ip . (I3.24p gets converted into a coupled semilinear system of real 
ODEs 

aA" + {{aw)'/w)A' - d($' - f3/a) 2 A + (7 + (3 2 /a)A + kA p+1 = 0, (3.26) 

a®" A + 2d($' - (5/ a) A! + {{an)' /m)A& - {{pw)' /w)A = 0, (3.27) 

which are seen to be the Euler-Lagrange equations SL/5A = and A~ 1 5L/5<& = of the 
Lagrangian (13.231) expressed in polar variables 

L = w(0 ( - a(0A' 2 - a{£)A 2 {& - fc)/a{£)) 2 + (t(0 + ^(0V«(0)^ 2 + 2k J A p+1 dA). 

(3.28) 

The phase rotation generator (I2.57P thereby becomes 

Y phas . = 9/9$, (3.29) 

producing a group of shift transformations $ — > $ + e with group parameter —00 < e < 00. 
Then the invariance prY p h as .£ = yields the multiplier equation 

wA(a<$>"A + 2d($' - p/a)A' + {{aw)' /vj)^' - {0w)'/w)A) = ^ {wA 2 {a& - p)), (3.30) 

which produces the first integral 

wA 2 {a& -P)=Ci (3.31) 

for the coupled semilinear system ( 13.26p ~l j3.27p . Through equation ( 13.3ip . $' can be elimi- 
nated in terms of A, so that the system reduces to a single real semilinear ODE 

aA" + {{aw)'/w)A + (7 + p 2 /a)A - {d 2 / W)yT 3 + kA p+1 = 0. (3.32) 

This ODE (I3.32p is the Euler-Lagrange equation 5L/5A = of the modified Lagrangian 

L=(L + 2(&-P{t)/a{t))C 1 ) 



-MO ( - «(£)- - 2 + (7(0 + ~P{tY/m)A 2 + {C^/a{0w{if)A- 2 + 2k I A p+1 dA) , 



<i) 2 )A- 2 + 2kJ A p+1 d 

(3.33) 

where C± is treated constant parameter. 

There are two cases, depending on C\ and p, for which the ODE (I3.32p can be solved. 
Firstly, if C\ — and p = —1 then the ODE becomes linear 

aA" + {{aw)'/w)A' + (7 + p 2 /a)A + k = (3.34) 

and its general solution will be given by special functions. Secondly, if C\ 7^ or p ^ — 1 then 
the ODE will reduce to quadratures by means of integrating factors arising from variational 
symmetries of the Lagrangian ( 13.331) . 
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Note that the set of solutions of the ODE (13.32)) for each value of G\ corresponds to the 
subset of solutions of the coupled semilinear system (13.26)) — ( T3.2T)) defined by the resulting 
level set of C\. Suppose a point symmetry generator 

Y = C{£)d/d£ + T{£)Ad/dA + u(£)d/d$ (3.35) 

admitted by the coupled semilinear system f l3.26p - fl3.27p leaves invariant the level set of Ci, so 
that prYCi = 2wA(T A-(A')(a® ' -/3)+waA 2 (cu-(&)' = when $' = (^+(Ci/ro)A- 2 )/a. 
Then, by reduction, the level-set ODE (13.321) will inherit a corresponding point symmetry 
generator 

Y = C(Od/d£ + T(g)Ad/dA. (3.36) 

Such symmetries are equivalent to conditional point symmetries of the form (I3.2p for the 
original complex 2nd order semilinear ODE (13. ip such that the generator preserves the first 
integral (13.311) given in terms of U, U by 

d =w(\\a(UU'- U'U) -PUU). (3.37) 

Any variational symmetries that are inherited in this way by the level-set ODE (13.32)) will 
yield the same integrating factors derived from the previous reduction method in section [37TI 
Now suppose that the level-set ODE (I3.32p admits a variational point symmetry that is 
not inherited from any point symmetry admitted by the coupled semilinear system (I3.26p - 
(I3.27p . Such a variational symmetry will be a hidden conditional symmetry which yields 
an additional integrating factor for the ODE (13. 32 p . The resulting reduction is most eas- 
ily carried out by utilizing canonical coordinates associated to the symmetry generator as 
follows. 

For a symmetry of the form (13.36)) . we change variables from (£, A) to canonical coordinates 
(z, B) given by 

Yz = 1, YB = 0, (3.38) 

where 

z = J(l/()d£, B = exp(- J(r/Od£)A. (3.39) 
The symmetry generator (13.36)) thereby becomes a ^-translation 

Y = d/dz (3.40) 
and thus the ODE (1 3 . 3 2 p is transformed into 

aB" + bB' + cB + qB' 3 + kB p+1 = (3.41) 

for B(z), where the coefficients a,b,c,q are real constants. Similarly, the Lagrangian (13.33)) 
becomes (modulo a total z-derivative) 

(L = L = exp{(b/a)z) ( - aB' 2 + cB 2 - qB~ 2 + 2k J B p+1 dB) , (3.42) 

dz 

where l/£ = — is the Jacobian of the transformation for the independent variable. Then 
Noether's theorem applied to the ^-translation generator (13.401) in characteristic form 



Y = -B'd/dB 
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(3.43) 



yields the multiplier equation 

prt(L) = - exp((6/a)z)^( - aB' 2 + cB 2 - qB~ 2 + 2k J B p+1 dB) 



B'^ + -^-(2(eM(b/a)z)aB')B'). 



(3.44) 



SB dz 

This equation (13.441) shows that b = is a necessary and sufficient condition for the z- 
translation generator (13.401) to be a variational symmetry of the ODE ( I3.32p . As an imme- 
diate consequence, the following result holds. 

Lemma 3. An ODE (13.411) has integrating factor B' iff b = 0. The resulting first integral 
is given by 

aB' 2 + cB 2 - qB- 2 + 2k J B p+1 dB = C 2 , (3.45) 

from which B(z) is determined by the quadrature 

dB 



±z + C 3 , (3.46) 



where 



H{B) = C 2 /a - (c/a)B 2 + {q/a)B- 2 - (2k /a) J B p+1 dB. (3.47) 

When 6 = 0, the quadrature (13.461) combined with the change of variable (13.391) will yield 
the general solution of the ODE (I3.32f) for A(C,). In this case the first integral ( 13.311) provides 
a quadrature for $(£), thereby yielding the general solution of the coupled semilinear system 
(I3.26p -( l3.27p . These quadratures for A(^) and $(£) can be written down in a fairly simple 
way. First, we find that the change of variable (I3.39P relating the Lagrangians (I3.42p and 
(I333|l gives 

A 2 IB 2 = exp (2 J (T/OO = (a(/ma) exp ((b/a) J (1/CR)- (3.48) 
Hence, from the quadrature (13.461) together with the variables (I3.39p . we see 

and 

A = {a(/w5t) 1/2 B (3.49b) 

for b = 0. Last, changing variables in the first integral ( I3.3ip and using the other first integral 
(I3.45p . we obtain 

$ = ±(d/o) J ^=^y + / G9/5)de + C A . (3.49c) 

The integrals (I3.49P yield the general solution of the coupled semilinear system (13. 26ft - 
(I3.27P for A(£) and $(^) under necessary and sufficient conditions for the reduced ODE 
( I3.32p to admit a variational symmetry (I3.36[) . This determines the corresponding solution 
£7(£) = A(£) exp(i$(^)) to the original complex semilinear Euler-Lagrange ODE (13. ip . ( I3.24p . 

Finally, suppose the level-set ODE ( I3.32p admits a hidden point symmetry that is not 
a variational symmetry. Such a hidden conditional symmetry can be used to obtain an 
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invariant solution of the ODE (I3.32p . which will determined a corresponding solution of the 
coupled semilinear system f !3.26|) — (13.2T|) through the first integral (13. 31 p as follows. 

Lemma 4. An ODE (I3.32p has a non-trivial invariant solution with respect to z -translations 
iffc^O or q ^ when p ^ —4, or c ^ and q ^ —k when p = —4. Then B' = yields 

= c + qB- A + kB p , (3.50) 

which is an algebraic equation determining B = const. ^ 0. 

Since B' = 0, relation (EQgj) yields 

A = (a(/wa) 1/2 Bexp ((b/2a) [ (1/CK). (3.51) 



Then, after changing variables from (£,A) to (z,B) in the first integral (I3.3ip . we find 
d$ ~ 

— — = (Ci/aB 2 ) exp(— (b/a)z) + ((3/a from which we obtain the quadrature 

$ = -CJibB 2 ) exp ((-b/a) J (1/C)de) + / + C 2 . (3.52) 

These expressions (I3.5ip - p.52p yield a solution of the coupled semilinear system ( 13. 26ft - 
( I3.27p . and hence U(£) = A(£) exp(i$(£)) gives a solution to the original complex semilinear 
Euler-Lagrange ODE ( 13. ip . ( I3.24p . corresponding to an invariant solution of the level-set 
ODE (I3.32p under a hidden conditional point symmetry ( I3.36[) . 

4. Solutions to the optimal translation-group ODE 

The reduction of order methods from section [3] will now be applied to the translation-group 
ODE ( 12 .5 Op arising by the reduction of the radial Schrodinger equation under its optimal 
subgroup of point symmetries ( I2.41bl) . 

As a preliminary step, we use polar variables U = Aexp(i$) to convert this [/(l)-invariant 
ODE (I2.50p into a semilinear system of real ODEs 

A" - A<$>' 2 + (n - ljr 1 ^' + vA + kA p+1 = 0, (4.1a) 

$" + 2A- 1 A'^' + (n - IjC 1 ®' = 0. (4.1b) 

The ODEs in this system (14. ip are the respective Euler-Lagrange equations SL/SA = and 
A^dL/SQ = of the {7(l)-invariant Lagrangian ( 12.521) expressed in polar variables, 

L = C~\-A 12 - A 2 & 2 + vA 2 + £-kA 2+p ). (4.2) 
In polar form, the point symmetries of ODE ( 12.501) listed in Proposition Q] consist of 

phase rotation Y p h as . = d/d&, (4.3a) 
scaling Y scal . = £d/di - (2/p)Ad/dA for v = 0, (4.3b) 
dilation Y dil . = £ 2 ~ n (£d/<9£ + (2 - n)Ad/dA) (4.3c) 
for z/ = 0, p = 2(3-n)/(n-2). 
Invariance of the polar Lagrangian (14. 2p under Y p h as . produces a first integral (13.311) . yielding 

$' = d^- n A- 2 . (4.4) 
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The polar system (14. ip thereby reduces to a single real semilinear ODE 

A" + {n- l)^- 1 ^ + vA + kA l+p - C\i 2 - 2n A- z = 0, (4.5) 

which is the Euler-Lagrange equation 5L/5A = of a modified Lagrangian (I3.33p . given by 

L = C~\-A' 2 + vA 2 + ^kA 2+p + Cli 2 - 2n A- 2 ). (4.6) 

Solutions of the ODE (14. 5 p for A(£) represent the level set C\ = const, of solutions 
(A(£), $(£)) to the polar system (14. ip . or equivalently the level set 

C x = \iC\UU' - U'U) = const. (4.7) 

of solutions U(£) to the translation-group ODE (I2.50p . 

Note that the level-set ODE (USD will be linear iff d = and p = -1. 

Proposition 2. In the nonlinear case C\ ^ or p ^ —1, the level-set ODE (14. 5p admits 
point symmetries only when 

v = 0. (4.8) 
For this case the admitted point symmetries consist of 

Y x =£d/d£ - (2/ p) Ad /OA, for C x = or p = 4/(n - 2), (4.9a) 

Y 2 =£ 2 -"(£/(n- 2)d/d£ -Ad/dA), for p = 2(3 - n)/(n - 2), (4.9b) 

Y 3 =2 ((C^/A;)^ 73 - £) 9/9^ - /or d ^ 0, p = -4, n = 4/3, (4.9c) 

Y 4 =2(Jfe - Cli 2 )id/di + {k- AC 2 ^ 2 ) Ad/dA, for C x ^ 0,p = -4, n = 0, (4.9d) 

Y 5 =^ 3 9/^ - 6(t 2 A - lQ/k)d/dA, for C 1 = 0,p=l,n= 16, (4.9e) 

Y 6 =9£ 2/3 d/d£ - (I2r 1/3 ^ + (4/A;)r 7/3 ) d/dA, for C 1 = 0,p=l,n= 13/3. (4.9f ) 

Since the point symmetries (j4.9aj) and (I4.9bp are the only ones also admitted by the polar 
system (14. ID . the remaining 4 point symmetries (I4.9cl) - (14.9fl) thus describe hidden conditional 
symmetries arising only for level-set solutions of the polar system (14.11) in the case v = 0, 
d ^ or p ^ -1. 

We now carry out the two reduction methods (cf sections 13.11 and 13. 2j) to obtain explicit 
solutions £/(£) for the translation-group ODE (I2.50p . The first reduction method can be 
applied only in the case v = 0, using the scaling and dilation symmetries (12.581) and (I2.59p . 
The second reduction method is applicable to the case v = by using the hidden conditional 
symmetries (I4.9cl) - (14.9fp . and the case v ^ by using the linearization which holds when 
C 1 = 0,p=-1. 

4.1. Scaling-symmetry quadratures for the translation-group ODE. The canonical 
coordinates of the scaling symmetry (12.581) are 

£ = ln£, V = i 2/p U (4.10) 

with v = 0. Hence the translation-group ODE (12.501) gets transformed into 

V" + (n - 2 - A/p)V + (2(2 - n)/p + A/p 2 )V + k\V\ p V = 0, (4.11) 

which is equivalent to the ODE (I2.60P describing static solutions of the radial Schrodinger 
equation. To apply the reduction Lemma [2] to ODE (14. lip , we note a = 1, b = n — 2 — 4/p, 
c = 2(2-n)/p+(2/p) 2 . 
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Then from Lemma [2](i) we obtain the invariant solution 

V = (2(n - 2 - 2/p)/kp) 1/p exp (r$), $ = const.. (4.12) 

This yields 

U = (2(n - 2 - 2/p)/A;p) 1/p r 2/p exp (r$), $ = const., (4.13) 

which is the scaling-invariant solution of ODE (I2.60p . 

Next, by Lemma [2In), we can obtain the general solution for V(z) in the case when 
n — 2 — A/p = under which the ODE (14.111) becomes 

V" - (2/p) 2 V + k\V\ p V = 0, p = 4/(n-2). (4.14) 

The quadratures (13.161) in polar variables V = Aexp(i<3>) are then given by 

<,A ± z + C 3 , $ = C 2 / - + C 4 , (4.15) 



(after renaming constants) where 

H(A) = Ci - C 2 2 yT 2 + (2/p) 2 A 2 - 2A; y A 1+p cL4. (4.16) 

These integrals (I4.15P cannot be evaluated generally to obtain explicit solutions for A(z) and 
$(z) in terms of elementary functions. Some special cases where explicit solutions can be 
derived are possible if we make a change of variables 

" (l/.)/^U y = A- (4.17) 



and, when C 2 ^ 0, 



v^p) 'J y/Hjy)' 



dA (W^SW y = A* (4.18) 



so that the respective expressions 

H x {y) = y^Hiy 1 '*) = C lV 2 ^ - C 2 2 y 2 ^ s + (2/pfy 2 - ~ky 2 ^ s (4.19) 

and 

H 2 {y) = y 2+2 / s H{y^ s ) = C iy 2+2 ^ s - C 2 2 y 2 + {2/pfy 2 ^ - ~ky 2+ ^ 5 for C 2 ± (4.20) 
are quadratic polynomials in y for some values of s and s, depending on p, C\, C 2 , where 

k = 2k/ (p + 2), p = A/{n - 2) ^ -2. (4.21) 
The required conditions from expression (I4.19p consist of 

2 + p/s = 0,1, 2; (4.22a) 

2-2/s = 0,l,2 ifC x ^0; (4.22b) 

2-4/s = 0,l,2 ifC 2 ^0. (4.22c) 

From expression (I4.20p . the required conditions are given by 

2 + 4/5 = 0, l,2and2 + (4 + p)/5 = 0,l,2 if C 2 ^ 0; (4.22d) 

2 + 2/5 = 0,1,2 if d ^ 0, C 2 ^ 0. (4.22e) 
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It is straightforward to solve these conditions for p, s, s. For each value found for p, only a 
single value each for s and s is needed to evaluate the integrals (I4.17P - (I4.18P . We thus find 

(4.23) 

2; (4.24) 
(4.25) 
(4.26) 

with n = 2 + A/p. By scaling and shifting y, we can then match the integrals (I4.17P - (14.181) 
to one of the forms 

/ ; = aTCCOS Hv), [ 7= = = arcsinh(y), (4.27) 

J y/y 2 - 1 J ^y 2 + 1 

{]>S In|y|, (4.28) 



Ci = and C 2 = : 


P 7^ 


-2,5 


= -P/2; 


Ci ^ and C 2 = : 


P = 


-M 


= 1; p = -4,S 


Ci = and C 2 ^ : 


P = 




= 2,s = -2; 


d ^ and C 2 + : 


P = 


-4,s 


= 2,s = -2; 



arcsin(y), (4.29) 



2/ 2 



(4 - 30) 

= y. (4.31) 



More generally, a similar method can be used to obtain additional explicit solutions for 
A(z) and <&(z) by requiring that the expressions fl4.19l) - fl4.20l) are either quartic polynomials 
in y, which will yield elliptic functions, or squares of quartic polynomials in y, which will 
yield elementary functions. (These solutions will be worked out elsewhere.) 

4.1.1. Quadrature for p ^ —2. From case ( 14.23p . we have s = —p/2, G\ = and C 2 = 0. 
The quadratures (13.161) are thus given by 

$ = C 4 (4.32) 

and 

± z + C 3 = {-2/p) [ dy y = A~ p l\ (4.33) 

J v / (4/P 2 )2/ 2 -2A;/(p + 2)' 

which can be matched to the form (|4.27|) . This yields 

A- p / 2 =p{2{p + 2)/k)' 1 / 2 cosh{±z + C 3 ), $ = C 4 , A;/(p + 2)>0, (4.34a) 

A- p/2 =p{-2{p + 2)/k)- 1/2 sinh{±z + C 3 ), $ = C 4 , k/{p + 2)<0. (4.34b) 

Hence we obtain two solutions of ODE (14. 14p for p ^ —2 

V = {2{p + 2)/{kp 2 )) 1 / p (cosh{±z + C 3 ))~ 2/p exp{iC 4 ), k/(p + 2) > 0, (4.35a) 

V = (-2(p + 2)/{kp 2 )) 1/p {smh{±z + C 3 ))- 2/p exp(iC 4 ), k/{p + 2) < 0. (4.35b) 

These solutions A4.35ap - fl4.35bp can be merged after we change variables ( I4.10p . giving the 
following result. 
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Proposition 3. The translation- group ODE (I2.60p has a solution 

U = (±8(p + 2)/(kp 2 )y/ p (C 3 e ± l/d)- 2 / p exp(id) (4.36) 

for 

P1 L -2,n = 2 + A/p,±k/(p + 2) > (4.37) 

with real constants C 3 , C4. 

4.1.2. Quadrature for p = — 1. From case (I4.24p . we have s — 1, C\ ^ and d = 0. The 
quadratures (13. 1 6[) are thus given by 

$ = d (4.38) 

and 

± z + C 3 = [ dy y = A, (4.39) 

J v/4(y - A;/4) 2 + d - P/4 V ; 

which can be matched to the forms (14.271) and (I4.28[) . This yields 

A = \ (k + (4Ci - £; 2 ) 1/2 sinh(2(±^ + C 3 ))) , $ = d, d > A; 2 /4, (4.40a) 

A = i (jfc + (k 2 - 4d) 1/2 cosh(2(±^ + C 3 ))) , $ = d, Ci < A; 2 /4, (4.40b) 

A = ±(A; + 4exp(2(±z + d))), $ = d, d = k 2 /A, (4.40c) 

A = i (/c — 4exp(2(±2 + C 3 ))) , $ = d, d = k 2 /4. (4.40d) 

Hence we obtain four solutions of ODE (14.141) for p = — 1 

V =\{k + (4Ci - A; 2 ) 1 / 2 sinh(2(±z + C 3 ))) exp(id), d > fc 2 /4, (4.41a) 

K = \ (k + [k 2 - 4Ci) 1/2 cosh(2(±^ + C 3 ))) exp(id), d < fc 2 /4, (4.41b) 

K = i(fc±4exp(2(2 + C 3 )))exp(iC 4 ), d = fc 2 /4, (4.41c) 

K = i(A;±4exp(2(-z + C 3 )))exp(id), d = ^ 2 /4. (4.41d) 

These solutions (I4.41al) - (14.41dp give the following result after we change variables (I4.10p . 

Proposition 4. The translation- group ODE (I2.60j) /ias solutions 

U = \ (ke + (±(d - A; 2 /4)) 1 / 2 (d^ 4 =F I/O,)) exp(id), d ^ A; 2 /4, (4.42) 

f/ = i(^ 2 + de 4 )exp(id), (4.43) 



/or 

rea/ constants C\, (7 3; d 



= i (*£ 2 + <%) exp(id) (4.44) 
p=-l,n = -2 (4.45) 
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4.1.3. Quadrature forp = —4. From cases (14.24fi — f)4.26p . we have s = 2, C\ ^ and G% = 0; 
s = 2, s = -2, and C 2 + 0. 

In the first subcase, the quadratures (I3.16P are given by 

$ = d (4.46) 

and 

±z + C 3 =f-= =^== = y = A 2 . (4.47) 
J y/(y + 2d) 2 + 4(fc - d ) 

Quadrature (14.471) can be matched to the forms (14.271) and (14. 28ft . This yields 

A 2 = -2d + 2(Jb - dY /2 sinh(±z + C 3 ), $ = d, C 2 < fc, (4.48a) 

A 2 = -2d + 2(C 2 - A;) 1/2 cosh(±z + C 3 ), $ = d, C 2 > fc, (4.48b) 

A 2 = -2d + exp(±z + d), $ = d, d 2 = ( 4 -48c) 

A 2 = -2d - exp(±z + d), $ = d, d 2 = k,d< 0. (4.48d) 

Hence we obtain four solutions of ODE ( 14. 14j) for p = —4 

V = (-2d + 2{k - d 2 ) 1/2 sinh(±^ + d)) 1/2 exp(id), d 2 < k, (4.49a) 

V = (-2d + 2(d 2 - A;) 1/2 cosh(±2 + C 3 )) 1/2 exp(id), d 2 > k, (4.49b) 

V = (-2d + exp{±z + d)) V2 exp(id), d 2 = k , ( 4 -49c) 

V = (-2d - exp(±z + d)) 1/2 exp(id), d 2 = k, d < 0. (4.49d) 

In the second subcase, the quadratures (13 . 161) are given by 

±Z + C 3 = f - dy — y = A 2 (4.50) 

and 

$ = C * - ° 2 \ m rJ\ J7 ^=T * = ^ ^ 4 - 51 ) 

J y/Mk - dV + 4dy + i 

Quadrature (I4.50p can be matched to the forms (I4.27[) and (I4.28p . which yields 

A 2 = -2d + 2(A;-d 2 -C , 2 2 ) 1/2 sinh(±2 + d), C 2 + C% < k, (4.52a) 

A 2 = -2d + 2(d 2 + d 2 - A;) 1/2 cosh(±z + d), d 2 + C l > k , (4.52b) 

A 2 = -2d + exp(±z + d), d 2 + Cl = k, (4.52c) 

A 2 = -2d - exp{±z + d), d 2 + C\ — k,d< 0. (4.52d) 
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Quadrature f)4.5ip can be matched to all of the forms fl4.27p - fl4.31l) . This yields 
$ =C 4 + \{k/C 2 2 - l)- 1/2 arcsinh ((Jfe - C\ - dT 1/2 (Ci + 2(k - C 2 2 )A- 2 )) , 

$ =C A - \{k/C 2 2 - l^^arccosh ((C 2 + C\ - k)- 1 ' 2 {C 1 + 2(k - C 2 2 )A- 2 )) , 

$ =d - |(1 - k/Cly 1 ' 2 arcsin ((C 2 + C? - A;)- 1 / 2 (C 1 - 2(C 2 - k)A~ 2 )) , 



(4.53a) 
(4.53b) 
(4.53c) 



$ = C 4 -(C 2 /(2C 1 ))ln|C , 1 A- 2 + l/2|, C 2 + C 2 = k, (4.53d) 
$ = d-(d/d)(d^T 2 + l/4) 1/2 , C 2 2 = fc, (4.53e) 
$ = d - C 2 A- 2 , C 2 = jfe, Ci = 0. (4.53f) 

Hence we obtain 7 more solutions of ODE f)4.14p for p = — 4 

K = (-2d + 2(Jfe - C\ - C 2 ) 1 ' 2 sinh(±^ + C 3 )) V2 

, /oWI . _2 , ,_ 1/2 . . / (fc - C 2 - C 2 ) 1 / 2 + d sinh(±z + C 3 ) \ \ 
x exp ^ - (l /2)(^ - 1) V 2arcsmh [ \ + l k _^_^ y ^ ±z + J z) ) ) , 

(4.54a) 

^ = (-2C 4 + 2(C 2 + C\ - A;) 1/2 cosh(±z + C 3 )) 1/2 

x exp ^iC 4 Ci/2K*/C 2 1) arccosh ^ _ (c| + ^ _ fc)1/2 + ^ j j , 

(4.54b) 

y = (-2d + 2(C 2 + C 2 - k) l/2 cosh(±z + d)) V2 

f r, , 1 /^2n-i/2 / (CI + C? " k ) 1/2 ~ Ci cosh(±z + d) \ \ 

x exp ^d - (1 /2)(1 - k/Cl) * arcsm ^ - ^ - ^ _ ^ - ^ j j , 

C*2 ^5 

(4.54c) 

V = (-2d + exp(±z + d)) 1/2 exp (iC 4 + i(d/(2d)) In |2d exp(=Fz - C 3 ) - l|) , 

Ci -I- C 2 = /c, 

(4.54d) 

V = (-2d - exp(±^ + d)) V2 exp (iC A + i(d/(2d)) In \2C X exp(=Fz - C 3 ) + l|) , 

-|- C 2 = fc, 

(4.54e) 

26 



V =(-2C, + 2| Cl | cosh( ±2 + C 3 )Y» exp U - .(^/(K?,)) ( ^ ^ + ^» ) " 2 ) 

V \-Ci + |Ci| cosh(±z + C 3 )/ / 

(4.54f) 

V = exp(i(±2 + C 3 )) exp (1C4 - iA; 1/2 exp(=F* - C 3 )) , C 2 = Jfe, d = 0. (4.54g) 

After we change variables (14.101) . these solutions (I4.49ap - fl4.49dp and (14 . 54al) ( [4 . 5 4g[ ) give 
the following result. 

Proposition 5. The translation- group ODE (I2.60P has solutions 

/ , ~ _ \ 1/2 

U = {-2C l i + {±{k-Cl)) 1 / 2 {C- i eTl/C,)) exp(itf 4 ), C\ ± k, (4.55) 

/ , \ 1/2 

U = ( ±2k 1/2 £ + C^ 2 ) exp(iC 4 ), k>0, (4.56) 

[/ = f±2A; 1 / 2 e + C 3 J exp(iC 4 ), k > 0, (4.57) 

[/ = (-2Ctf + (A; - C 2 - C 2 ) 1 ' 2 ^ 2 - \/C*f ^ 



x exp U - H k,d - l)-V 2 arc S1 nh( 2 <* = 3 = ^ + ^ = ^ ) V 
V 2 V-2C 1 e + (A;-C , 1 2 -C , 2 2 ) 1 /2(C3e-l/C3)V 

(4.58) 

U = (-2Ctf + (C 2 + C\ - k) 1/2 (c 3 e + i/c* 3 )) 1/2 



x exp U - W ,Cl - l)-V 2 arccosh( 2 ^ 2 + ^ 2 - fc ) 1/2 ^^^ 2 + 1 ^ 
V 2 V2C 1 e-(C 2 + C 1 2 -A;)V2(C3e + l/C3) 



(4.59) 

[/ = (-2Ctf + (C 2 2 + C 2 - fc) 1 / 2 (^ 2 + l/Ca)) 172 

x exp 1C4 — — fc/Co 7 arcsm ■ — ~ — , 

V 2 \-2C^ + (C 2 + C 2 -ky/^C 3 e + l/C 3 )'>J 

(4.60) 

U = (-2&Z + C 3 e 2 ) V2 exp (i<7 4 + i(C a /(2Ci)) In ^C^" 1 - C'al) , C 2 + C 2 = fc, 

(4.61) 

U = (-2C& + <7 3 ) V2 exp (1C4 + i(C 2 /(2C7i)) In |2d£ - C 3 |) , C 2 + C 2 = k, (4.62) 
17 =(C x /C*yi\C& - 1) exp (iC 4 - i^/d) (^37)) . * > °' ( 4 - 63 ) 

[/ =C* 3 1/2 £ exp (1C4 - i^/OOr , fc > 0, (4.64) 
U =Cl' 2 exp (1C4 - i(A; 1/2 /(7 3 )^ , fc > 0, (4.65) 
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for 

p= -4,n = l (4.66) 

with real constants C\, C 2 , C3, C±, C* 4 . 

4.2. Dilation-symmetry quadratures for the translation-group ODE. The canonical 
coordinates of the dilation symmetry (12.59)) are 

z = C~ 2 , V = C~ 2 U (4.67) 

with v — 0. In this case the translation-group ODE (I2.60P gets transformed into 

V" + k(l+p/2) 2 \V\ p V = 0, p = 2(3 -n)/(n- 2), n ^2, 3. (4.68) 

To apply Lemma |2]to this ODE pg| . we note a = 1/(1 + p/2) 2 , 6 = 0, c = 0. 

From Lemma [2(i) we see that the invariant solution is trivial, since the only root of F is 
\V\ =0. Thus ODE (12.601) has no non-trivial dilation-invariant solution. 

Next, we can obtain the general solution for Viz) from the quadratures (13.161) . since 
the conditions of Lemma [2](ii) are directly satisfied for the ODE (14.681) . In polar variables 
V = v4exp(i$), these quadratures are given by 

lJA ± z + C 3 , $ = C 2 [ —^== + C 4 , (4.69) 



where 



= Ci - Cf - 2fc(l + p/2) 2 J A 1+p dA. (4.70) 

These integrals (14.69)) cannot be evaluated generally to obtain explicit solutions for A(z) and 
<&{z) in terms of elementary functions. Some special cases where explicit solutions can be 
derived are possible if we change variables (I4.17I) - (I4.18I) to make the expressions 

H x {y) = y 2 - 2 / s H(yV*) = C iy ^ s - C 2 2 y 2 ~^ - % 2+ ^, (4.71) 

H 2 (y) = y 2+2 ^H(y 1 ^) = C lV 2+2/rs - C 2 y 2 - ky 2+ ^ +p)/S for C 2 ? (4.72) 

into quadratic polynomials in y for some values of s and 5, depending on p, Ci,C 2 , where 

k = 2k/(p + 2), p = 2(3 -n)/(n- 2) ^-2. (4.73) 

The required conditions from expressions (14.711) and (14.721) are respectively given by (I4.22a[) . 
(l4T22bl) . (I4722c1) . (l4722e]) and 

2 + (4 +p)/5 = 0,1, 2 ifC 2 ^0. (4.74) 
Solving these conditions for p, s, s in the same way as for the previous case, we find 



Ci = and C 2 = : 


P 7^ 


-2,s 


= -P/2; 


(4.75) 


d ^ and C 2 = : 


P = 


-1,5 


= 1; p=-4, S = 2; 


(4.76) 


d = and C 2 ^ : 


P = 


-4,s 


= 2, 5=-l; p = -8, s = 4, 5 = 2; 


(4.77) 


d ^ and C 2 ^ : 


P = 


-4,5 


= 2,5= -1; 


(4.78) 



with n = 2(p + 3)/(p + 2). 

We can then match the integrals (I4.17l) -( l4.18l) to one of the forms (14.271) ( 14.311) after 
scaling and shifting y. 
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As noted previously, additional explicit solutions for A(z) and $(-2) can be derived by 
requiring that the expressions ( 14 . 7 1 [) — ( 14 . 72 j) are either quartic polynomials in y, which will 
yield elliptic functions, or squares of quartic polynomials in y, which will yield elementary 
functions. (These solutions will be worked out elsewhere.) 

4.2.1. Quadrature for p 7^ —2. From case (14.75 p . we have s = —p/2, C\ — and C2 = 0. 
The quadratures ( 13.1 6p are thus given by 

$ = C 4 (4.79) 

and 

±a + C ^- (2/rt / v-% + 2)/2 - » = ^ (4 - 80) 
which can be matched to the form (14.3 ip . This yields 

A- p ' 2 = {-kp 2 {p + 2)/8) 1/2 {±z + C 3 ), $ = C 4 , k{p + 2)<0. (4.81) 

Hence we obtain a solution of ODE ( I4.68P for p ^ —2 

V= {-kp 2 {p + 2)/8)- 1/p {±z + C 3 )- 2/p exp{iC±), k{p + 2)<0. (4.82) 

This solution ( 14 . 8 2 [) gives the following result after we change variables (I4.67p . 

Proposition 6. The translation- group ODE ( I2.60P has a solution 

U = (-kp 2 {p + 2)/8)- 1/p {£ + C , 3^ n )~ 2/p exp(iC 4 ), k{p + 2) < (4.83) 

for 

p^ -2,n = 2(p + 3)/(p + 2) (4.84) 

with real constants C3, C4. 

4.2.2. Quadrature for p = —1. From case (14. 76 p . we have s — 1, C\ ^ and C2 = 0. The 
quadratures ( 13 . 1 6[) are thus given by 

$ = C 4 (4.85) 

and 

± z + C 3 = [ ; dy y = A, (4.86) 

J ^C x -(kj2)y K } 

which can be matched to the form (I4.30|) . This yields 

A = {2C x /k) - (k/8)(±z + C 3 ) 2 , $ = C 4 . (4.87) 
Hence we obtain a solution of ODE ( 14 .681) for p = —1 

V = ({2C x lk) - (k/8)(±z + C 3 ) 2 ) exp(iC 4 ). (4.88) 
This solution ( 14 . 8 8 1) gives the following result after we change variables (14.671) . 
Proposition 7. The translation- group ODE ( I2.60p has a solution 

U = ((2d/fc)r 2 - (*/8)(£ + C 3 /0 2 ) exp(iC 4 ) (4.89) 

for 

p=-l,n = 4 (4.90) 

with real constants C\, C 3 , C 4 . 
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4.2.3. Quadrature for p = —4. From (14. 76 1) - (14. 78 p . we have s = 2, s = — 1 and C 2 ^ 0. 
The quadratures ( 13 . 1 6[) are thus given by 



± z + d = (1/2) / - » y = A , (4.91) 
J y/dy + k-C 2 

which can be matched to the forms fl4.30[) — (14.31 j) . and 

$ = C 4 - C 2 f dV =, y = A-\ (4.92) 

which can be matched to all of the forms ( I4.27l) -( ]4.29l) . This yields 

^ — (Cf - k)/Ci + Ci(±z + C 3 ) 2 , d^O, (4.93a) 

A 2 = 2(fc - C|) 1/2 (±^ + C 3 ), C| < k,Ci = 0, (4.93b) 

and 

$ =d - C 2 (Jfe - C 2 2 )- 1/2 arcsinh (((jfe - Cf yCi) 1 ^ -1 ) , C 2 2 < k,C 1 > 0, (4.94a) 

$ =d - C 2 [C\ - A;)~ 1/2 arcsin (((C 2 2 - jfeVd) 172 ^" 1 ) , C 2 > jfe, d > 0, (4.94b) 

$ =C 4 - C 2 (ife - C 2 )- 1/2 arccosh (((C 2 - k) j Cx)^ 7, A' 1 ') , C 2 < k, d < 0, (4.94c) 

$ =d - (d/Cj 72 )^ 1 , C 2 = jfe, d > 0, (4.94d) 

$ =C 4 -C 2 (A;-C 2 2 )- 1/2 ln|A- 1 |, C 2 < k,C x — 0. (4.94e) 
Hence we obtain 5 more solutions of ODE (14. 14h for p = —4 



V = {{Cl-k)/C l + C l {±z + C,) 2 ) 1 ' 2 

x exp fid - iC 2 {k - C 2 )- 1/2 arcsinh(^(-l + C 2 (±z + C 3 ) 2 /(k - C 2 )) 



1/2 

CI < jfe, d > 0, 

2\l/2 



(4.95a) 

V = ((C 2 - £;)/d + C x {±z + C 3 ) 2 ) ] 

x exp fid - iC 2 (C£ - ky 1 ' 2 arcsin ((1 + C 2 (±z + d)7(<?| - £0)~ 1/2 )) , (4.95b) 

C 2 >k,d> 0, 

V = ((C 2 2 -k)/C 1 + C 1 (±z + C 3 ) 2 ) 1/2 

x exp ^d - iC 2 {k - C 2 )~ 1/2 arccosh((l + C?(±* + C 3 ) 2 /{C% - £0)~ 1/2 )) , (4.95c) 



C\ < k, Ci < 0, 

tl/2i 



V =Ci /2 | ± 2 + d| exp (1C4 =F iC/c^Vd)! ± * + dl -1 ) , C 2 2 = k, d > 0, (4.95d) 

V =(4(Jb - C 2 2 )) l/ \±z + C 3 ) 1/2 exp (id + \iC 2 {k - C 2 )- 1/2 In | ± z + d 



C 2 < jfe, d = 0. 



(4.95e) 



7 2 

After we change variables (14. 6 7ft , these solutions (14.95ajl — (I4.95e[) give the following result. 
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Proposition 8. The translation- group ODE (I2.60P has solutions 

1/2 



U = ((CQ-k)e/C l + C 1 {l + C s £f 

x exp (1C4 - iC 2 {k - Cf )- 1 / 2 axcsinh((-l + C 2 ^ 1 + C 3 ) 2 /(k - C 2 ))- 1 ^ , ( 4 - 96 ) 



C 2 <k,d> 0, 

1/2 



^ = ((^ 2 2 -^ 2 /Ci + Ci(i + c'30 2 ' 

x exp ^iC 4 - iC 2 (C 2 2 - k)- 1 ' 2 arcsin ((1 + C?^" 1 + C 3 ) 2 /(C% - fc))~ 1/2 )) , ( 4 - 97 ) 
Cf > Jfc, Ci > 0, 

/ \ 1/2 

u = [(c 2 2 -k)e/c 1 + c 1 (i + c 3 o 2 ) 

x exp ^iC 4 - iC 2 (k - C 2 )- 1/2 arccosh((l + C?^" 1 + C 3 ) 2 /(C 2 2 - fe)) -1 / 2 ^ , ( 4 ' 98 ) 



Cf < fc, d < 0, 



17 =(4(* - (%)YI\Che ± 1/2 ^p (i<7 4 + \iC 2 {k - Cj)- 1 ' 2 In |C 3 ± C'l) , C 2 2 < k, 

(4.99) 

U =d /2 \l + exp (iC 4 =f i^/COir 1 + C'sr 1 ) , C a > 0, (4.100) 

for 

p= -4,n=l (4.101) 

with real constants C\, C 2 , C 3 , C 4; C 4 . 

4.2.4. Quadrature forp = —8. From case (I4.77p . we have s = 4, 5 = 2, Ci = and C 2 7^ 0. 
The quadratures (13 . 1 6[) are thus given by 

±z+c - (1/4, /^ra- ^ 4 (4102) 



and 



Quadrature (14.1021) can be matched to the form (14.301) . which yields 

A 4 = 3k/Ci - 4Ci{±z + C 3 ) 2 , k > 0. (4.104) 
Quadrature (14.1031) can be matched to the form (14.291) . This yields 

$ = C 4 + (1/2) arcsin ({3k)^ 2 C 2 A 2 ) , k > 0. (4.105) 
Hence we obtain a solution of ODE (I4.68P for p = —8 

V = {Zk/Cl - ACl{±z + C 3 ) 2 ) 1/4 exp (1C4 + ±iarcsin((l - 4C 2 4 (±^ + C 3 ) 2 /(3A;)) 1/2 )) , 
k > 0. (4.106) 
This solution (14. 106ft gives the following result after we change variables (14. 67ft . 
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Proposition 9. The translation- group ODE (I2.60p has a solution 

1/4 



u = ((3k/c 2 2 )e /3 -ic 2 2 (e /3 +c 3 z 



2/3 \ 2 



f , \\ (4 - 107) 

x exp 1 1C4 + |iarcsin ((1 - 4C 2 4 (r V3 + C3) 2 / (3k))^ 2 j J, fc > 

for 

p=-8,n = 5/3 (4.108) 

rea/ constants C2, C3, C4. 

4.3. Conditional-symmetry quadratures for the translation-group ODE. From 
Proposition [21 the hidden conditional symmetries ( I4.9cp - (14.9fj) with v — can be natu- 
rally split up into two types: symmetries (I4.9cl) . (14.9dl) hold for p = —4, C\ 7^ and have 
the form (13.361) to which Lemmas [3] and [4] can be applied; symmetries (I4.9el) . (14.9f|) hold for 
p = 1, C\ = and have a slightly more general form 

Y = ai)d/di + (Tx(£)A + T {Z))d/dA, (4.109) 

for which the methods that underlie Lemmas [3] and H] still can be used. 

4.3.1. Quadratures for p = —4, C\ 7^ 0. Point symmetry (|4.9c[) of the level-set ODE (I4.5[) 
has the canonical coordinates 

z = (-3/4) HC^/k - £ 2/3 ), 5 = [C x 2 jk - ^ 3 y 3 ^A (4.110) 

with n = 4/3, z/ = 0. Hence the Lagrangian ( 14.61) for ODE ( 14.51) gets transformed into 

L = exp(-(2/3)z)(-(l/2) J B /2 - (l/6)5 2 + 2kB~ 2 ), (4.111) 

yielding the transformed ODE 

\B" - \B' - \B - 2kB~ 3 = 0. (4.112) 

Since the coefficient of B' is non-zero, the reduction shown in Lemma [3] cannot be applied 
to ODE (14. 112)) . Instead, from Lemma HI an invariant solution can be derived. This yields 

B = (-12fc) 1/4 . (4.113) 

Substituting this solution into expressions ( 13 .5 If) ( 13.52)) . we obtain 

A = (-Uk^C^/k - £ 2/3 ) 3/4 , $ = (V3d/2)(k^ 3 - d 2 )- 1/2 + C 2 , (4.114) 

which is the corresponding invariant solution of the coupled semilinear system (13.26)) — (13.27)) 
with p = —4, n = 4/3, v = 0. 

Point symmetry ( 14.9dj) of the level-set ODE ( 14. 5 p has the canonical coordinates 

z = (1/(4*)) ln(e 2 /(Ci 2 f - *0), B = T 1/2 m 2 - k)~^A (4.115) 
with n — 0, v — 0. Hence the Lagrangian ( 14. 6 j) for ODE (14.51) gets transformed into 

L = exp{-2kz){{l/2)B' 2 + {3k 2 /2)B 2 - 2B' 2 ). (4.116) 
This yields the transformed ODE 

- \B" + kB' + \k 2 B + 2B- 3 = 0, (4.117) 
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whose coefficient of B' is non-zero. Again, while the reduction shown in Lemma [3] cannot be 
applied, an invariant solution to ODE (I4.117P can be derived from Lemma HJ yielding 

B = (-3A; 2 /4)" 1/4 . (4.118) 
Substituting this solution into expressions ( I3.5ip -( 13.52p . we obtain 

A = (3A; 2 /4)- 1/4 e 1/2 (A; - Ci 2 £ 2 ) 3/4 , $ = -(v / 3C 1 /2)(A;r 2 - Ci 2 )~ 1/2 + C 2 , (4.119) 

which is the corresponding invariant solution of the coupled semilinear system ( 13.26l) -( |3.27l) 
with p = —4, n = 0, v = 0. 

These two solutions ( I4.114p and (I4.119P give the following result. 

Proposition 10. The translation- group ODE ( I2.60p has solutions 

U ={-Ylk) x l\C x lk - e 2/3 ) 3/4 exp (i(V3/2)((k/C 1 )£ 2/3 - 1)- 1/2 + iC 2 ) (4.120) 

for 

p= -4,n = 4/3, k <0 (4.121) 

and 

U =(4A;/3) 1/4 e 2 (r 2 - C'i/A;) 3/4 exp ( - 1(^/2) ((k/C^C 2 - 1)" 1/2 + iC 2 ) (4.122) 

for 

p=-4,n = (4.123) 

with real constants C\, C<i- 

4.3.2. Quadratures forp = 1, C\ = 0. Point symmetry (I4.9ej) of the level-set ODE (14.51) has 
the canonical coordinates 

z = -(i/2)r 2 , B = eA-{2A/k)e (4.124) 

with n — 16, Cx — 0, v — 0. Hence the ODE (14.51) gets transformed into 

B" + kB 2 = 0. (4.125) 

This ODE (14.1251) admits the integrating factor B', leading to the quadrature 

dB 



, = = ±z + C 3 , (4.126) 

ySC 2 - (2k/3)B 3 

which is an elliptic integral when Ci ^ 0. (This solution will be worked out in detail 
elsewhere.) In the case C*2 = 0, the integral ( 14.1261) can be evaluated explicitly, yielding 

B = (-6/k)(±z + C 3 y 2 . (4.127) 

Hence, from the canonical coordinates (14.1241) and the first integral (13 . 3 1 [) . we obtain 

A=(9QC 3 /k)(C 3 e±l)/(2C 3 e±l) 2 , $ = C 4 , (4.128) 

which is the solution of the coupled semilinear system ( 13.26j) (13. 271) with p = 1, n = 16, 
v = 0, and C\ — C-2 — 0. Note that the invariant solution of the ODE (14. 1251) . given by 
B' = 0, is simply B = 0. This solution also arises as the special limit C 3 — > oo in the solution 
( I4.127p . The corresponding invariant solution of the coupled semilinear system (13. 26j) — f T3.27f) 
is given by 

A = 24/(^ 2 ), $ = C 4 , (4.129) 
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which is non-trivial due to the shift term in the change of variables (I4.124p for B in terms 
of A. 

Point symmetry ( 14.9fj) of the level-set ODE (14.51) has the canonical coordinates 

z = (1/3)£ 1/3 , B = e 4/3 A - (2/(3£;))r 2/3 (4.130) 
with n = 13/3, v = 0, G\ = 0. Hence the ODE (14.51) gets transformed into 

±B" + kB 2 = 0. (4.131) 
This ODE (14.1311) admits the integrating factor B', leading to the quadrature 

dB 



= = ±z + C 3 . (4.132) 

V81C 2 - 54fc£ a V 1 

When C 2 7^ 0, this quadrature is an elliptic integral. (The resulting solution will be worked 
out in detail elsewhere.) In the case C 2 = 0, the integral ( 14.132)) can be evaluated explicitly, 
yielding 

B = -(2/27k)(±z + C 3 )- 2 . (4.133) 
The canonical coordinates (14.124)) and the first integral (13.311) then give 

A = (2C 3 A)r 2 (3C 3 ± 2£ 1/3 )/(3C 3 ± e /3 ) 2 , $ = C 4 , (4.134) 

which is the solution of the coupled semilinear system (13.26)) (13.27^ with p — 1, n = 13/3, 
v = 0, and C\ — C 2 = 0. Note that the invariant solution of the ODE ( 14 .131)) is again B = 
which also arises as the special limit C 3 — > oo in the solution (14. 133)) . The corresponding 
invariant solution of the coupled semilinear system ( 13.26)) — ( ]3 .27)) is given by 

A = 2/{3k£ 2 ), $ = C 4 . (4.135) 

The solutions ( 14TT281) . (14129]) . (Q3"i|) . (I4TT35]) give the following result. 

Proposition 11. The translation- group ODE ( 12.60)) has solutions 

U =(24/A;)r 2 exp(iC , 4 ), (4.136) 



for 
and 



U =(96A)(£ 2 + C 3 )(2e + C 3 r 2 exp(iC7 4 ), (4.137) 

p= l,n = 16 (4.138) 

C/ =(2/(3A;))r 2 exp(iC 4 ), (4.139) 

U =(2/k)C 2 (2C 3 e /3 + S)(C 3 e /3 + 3)- 2 exp(iC 4 ), (4.140) 



for 



p = l,n =13/3 (4.141) 

TOi/i rea/ constants C 3 , C 4 . 
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4.4. Linearization of the translation-group ODE. When C\ = and p = — 1, the 

level-set ODE (14.51) for becomes linear, 

A" + (n-l)C 1 A' + vA + k = 0. (4.142) 

The general solution of this ODE splits into three cases. 

For the case v > 0, we can use a change of variable A = £ 1_n / 2 .B, z = which 
transforms the ODE (14.1421) to an inhomogeneous form of Bessel's equation 

z 2 B" + zB' + {z 2 - (|n - lf)B + ku- (2+n ^ 4 z 1+n/2 = 0. (4.143) 

This equation has the general solution 

B = C 2 J\ n „ 2 \ /2 {z) + C 3 Y ]n _ 2]/2 (z) - kiy-^^z^ 2 , (4.144) 

where J M (^) and Y^z) are the Bessel functions of the first and second kinds, respectively. 
Hence the general solution of ODE ( 14.1421) for A(^) with v ^ is given by 

A = e~ n/2 (C 2 J|„_ 2 |/2(V^0 + C 3 l|n-2|/2(V^0) " k/v. (4.145) 

For the case v < 0, we can use a change of variable A = £ 1 ~ n / 2 B, z = a/— z/£ which transforms 
the ODE (I4.142p to an inhomogeneous form of a modified Bessel equation 

2 2 5" + 25' - (z 2 + (|n - l) 2 )B + A;(-z/)-( 2+n )/S 1+n/2 = 0. (4.146) 

This equation has the general solution 

B = C 2 /|n-2|/ 2 (z) + C 3 K ln . 2l/2 (z) + k(-up 2+n ^z\ n - 2 \/ 2 , (4.147) 

where I^{z) and K^(z) are the modified Bessel functions of the first and second kinds, 
respectively. Hence the general solution of ODE (14. 142j) for A(^) with v ^ is given by 

A = e~ n/2 (C 2 /| n _ 2 | /2 ( v / ^0 + C 3 K ln _ 2l/2 (V^O) - k/v. (4.148) 
Since C\ = 0, the first integral (13.311) for $(£) reduces to $' = 0, which yields 

$ = const.. (4.149) 
For the case v = 0, the ODE (14.1421) can be integrated directly to get 

A= -{k/2n)i 2 + C 3 i 2 ~ n + C 2 , n^0,2, (4.150) 

A = -(k/4)£ 2 + C 3 \n£ + C 2 , n = 2, (4.151) 

A= (-k/2)t 2 \n£ + C 3 Z 2 + C 2 , n = 0. (4.152) 
From the first integral (I3.31f) for $(£), we again have 

$ = const.. (4.153) 

Hence we obtain the following result. 
Proposition 12. The translation- group ODE (I2.50p has solutions 

U = (^- n/2 (C 2 J| n _ 2 , /2 (^) + C 3 ^_ 2 | /2 (V^0) - k/u) exp(iC 4 ) (4.154) 

for 

p = -l,z/ > 0,Ci = (4.155) 

and 

C/ = (e 1_n/2 (C 2 /|n- 2 |/ 2 (v / ^0 + O^n-siM^O) - k/u) exp(iC 4 ) (4.156) 
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for 




p = —i-jV < U,Oi = U 


[4:. LQl ) 


and 










TT - 
U - 


-[ — {K/Zn)l; + o 3 £ + o 2 J exp(io 4 J 


(A 1 Cifi^ 
yt. loo ) 


for 




P = — 1,B f U,flf = U,Oi = U 


//i i £r\\ 


and 










TT 

u 


— \~ l fc / 4 A + o 3 m^ + o 2 J exp(io 4 j 


(^4. lOU J 


for 




p — — 1, 71 — i/ — U,Oi — u 


(^4. lOlJ 


and 










U 


= ((-k/2)e\n£ + C 3 e + C 2 )exp(iC 4 ) 


(4.162) 


for 




p = = 0,i/ = 0,Ci = 


(4.163) 


with real constants Ci, 


c 3 , 


C 4 . 





5. Solutions to the optimal scaling-group ODE 

Next the reduction of order methods from section [3] will be applied to the scaling-group 
ODE (12.461) arising by the reduction of the radial Schrodinger equation under its optimal 
subgroup of point symmetries (12.41 cp . 

To proceed, we use polar variables U = Aexp(i$) to convert this [/(l)-invariant ODE 
(I2.46P into a semilinear system of real ODEs 

Ai 2 A" - Ai 2 A& 2 - 2(n - 4 - 4/p)£A' + (1 + 4n£)A& (5.1a) 

- (/i 2 - (4 - 2n) Ip - A/p 2 )A + kA 1+p = 0, 

4fV + g^A^A'®' - (1 + 4/i£)A- 1 A / - 2(n - 4 - 4/p)£$' + //(n - 2 - 4/p) = 0. (5.1b) 

As shown by the results stated in Proposition [T], the point symmetries of this system (15.11) 
consist of only phase rotations 

Y phas . = d/d<$>. (5.2) 

In the case of the conformal power p — 4/n, the ODEs in the polar system 05. ip are the 
respective Euler-Lagrange equations SL/SA = and A~ 1 5L/5§ = of the [/(l)-invariant 
Lagrangian (12.471) expressed in polar variables, 

L = _4£2 ( ^ + A 2 ( <I>,_ ( i^ p = 4/n. 

(5.3) 

Invariance of the polar Lagrangian (I5.3P under Y p h as . produces a first integral (13.311) . yielding 

= |r 2 (C , iA- 2 + l + 4/iO, P = 4/n. (5.4) 
The polar system 05. 1 p thereby reduces to a single real semilinear ODE 

4£ 2 A" + 8£A' - \C 2 C 2 A- 3 + ( ^r 2 + l/^r 1 + 4^i)A + kA 1 ** = 0, p = 4/n, (5.5) 
which is the Euler-Lagrange equation 5L/5A = of a modified Lagrangian (I3.33p . given by 

L = _ A ? A ' 2 + \C\C 2 A- 2 + (^r 2 + + ^)A 2 + ^r 2 kA^ 2 , p = 4/n. (5.6) 
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Solutions of the ODE (15. 5p for A(£) represent the level set C\ = const, of solutions 
(A(£), $(£)) to the polar system (15. ip . or equivalently the level set 

d = 2i£ 2 {UU' - U'U) - (§ + 2^)\U\ 2 = const. (5.7) 

of solutions U(£) to the scaling-group ODE (12.461) . 

Note that the level-set ODE (15. 5p will be linear iff G\ = and p = — 1. 

Proposition 13. In the nonlinear case C\ ^ or p ^ —1, the level-set ODE (15. 5ft admits 
no point symmetries. 

Since both the level-set ODE (15. 5ft and the polar system (15.11) have no point symmetries, 
only the second reduction method (cf section I3.2j) is applicable by using the linearization of 
ODE (15.51) which holds when G\ — 0, p — — 1. 



5.1. Linearization of the scaling-group ODE. In the case G\ = andp = — 1 (n = —4), 
the level-set ODE (15.51) for A(£) becomes linear, 

A? A" + 8£A' + (^r 2 + §^ _1 - 8) A + k = 0. (5.8) 
The general solution of this ODE is given by 

A =C 2 M^ /2m (\iC l ) + C 3 M_ I/l/2 ,_3/ 2 (i^- 1 ) 

+ (ifc/3)(M_ W2i _ 3 / 2 (|ir 1 ) f M_ 1# , As/a (^- 1 ) de (5.9) 

- M-^/aCK" 1 ) / M^fr-wijpC 1 ) d£), 

where M\ v (z) is the Whittaker function. 

Since G\ = 0, the first integral (13.311) for $(£) reduces to $' = |£~ 2 + §/i£~\ which yields 

$ = -|r 1 + |/ilne + C , 4- (5.10) 
Hence we obtain the following result. 
Proposition 14. The scaling-group ODE (I2.46P has a solution 

U =(c 2 M_ W2 , 3/2 (iir 1 ) + CaM-^A-aMK- 1 ) 

+ (ifc/3)(M_ W2 ,_ 3 / 2 (iir 1 ) y M_ 1/i/2j3 / 2 (jir 1 )^ 

- M_ W2i3/2 (iir 1 ) / M_ W2) _ 3/2 (iie- 1 ) dCj) exp(-|ir 1 + |i/ilne + iC 4 ) 

(5.11) 

for 

p=-l,n = -4,Ci = (5.12) 

witt rea/ constants C 2 , C3, C4. 
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6. Solutions to the optimal conformal-group ODE 

The reduction of order methods from section [3] will be considered last for the conformal- 
group ODE ( I2.55P arising by the reduction of the radial Schrodinger equation under its 
optimal subgroup of point symmetries (I2.41dj) . 

We begin by using polar variables U = v4exp(i$) to convert this [/(l)-invariant ODE 
( I2.55P into a semilinear system of real ODEs 

Ai 2 A" - 4£ 2 A$' 2 + 8£A' + (kC 1 - \C 2 + n(l - \n))A + kA l+4/n = 0, (6.1a) 
$" + 2A~ 1 A'& + = 0. (6.1b) 

The ODEs in this system (16.11) are the respective Euler-Lagrange equations SL/SA = and 
A^SL/SQ = of the £/(l)-invariant Lagrangian (12.561) expressed in polar variables, 

L = -A? A' 2 - 4£ 2 A 2 $' 2 + (kC 1 - |r 2 + n{l - \n)) A 2 + -^kA 2+4 / n . (6.2) 

As shown by the results stated in Proposition [H the point symmetries of the polar system 
(16. ip consist of only phase rotations 

Y phas . = d/d§. (6.3) 

Invariance of the polar Lagrangian (I6.2p under Y phas produces a first integral (13.311) . yielding 

$' = \dC 2 A- 2 . (6.4) 

The polar system (16. ip thereby reduces to a single real semilinear ODE 

4£ 2 A" + 8£A' - \ClC 2 A- z + (kC 1 - \C 2 + <l - \n)) A + kA 1+A/n = 0, (6.5) 

which is the Euler-Lagrange equation SL/5A = of a modified Lagrangian (I3.33p . given by 

I = ^a 12 + \C\C 2 A- 2 + (kC 1 - \C 2 + n(l - in)) A 2 + ^- 2 kA 2+i ' n . (6.6) 

Solutions of the ODE (16. 5p for represent the level set C\ = const, of solutions 

(A(£), $(0) to the polar system (16. ip . or equivalently the level set 

C x = 2iC 2 (UU' - U'U) = const. (6.7) 

of solutions U(£) to the conformal-group ODE (I2.55p . 

Note that the level-set ODE Qg2J will be linear iff C x = and p = -1 (n = -4). 

Proposition 15. In the nonlinear case C\ ^ or p ^ — 1 (n ^ — A), the level-set ODE 
(16. 5p admits no point symmetries. 

Since both the level-set ODE (16. 5p and the polar system (16.11) have no point symmetries, 
only the second reduction method (cf section I3.2p is applicable by using the linearization of 
ODE K52J) which holds when d = 0, p = -1 (n = -4). 

6.1. Linearization of the conformal-group ODE. In the case C\ = and p = —1 

(n = —4), the level-set ODE (16.51) for A(£) becomes linear, 

Ai 2 A" + 8£A' + « _1 - lr 2 - 8) A + k = 0. (6.8) 
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The general solution of this ODE is given by 

A =C 2 M K/2i3/2 (ir 1 ) + OiAk/a.-aMir 1 ) 

+ (fc/6)(M K/2i3/2 (ir 1 ) J M K/2 ^ /2 {\C l )d^-M K/ ^ /2 {\C l ) j M^a/adr 1 )^), 

where M^ ^z) is the Whittaker function. 

Since C\ — 0, the first integral (I3.3ip for $(£) reduces to $' = 0, which yields 

$ = const.. (6.10) 

Hence we obtain the following result. 
Proposition 16. The conformal- group ODE (I2.55P has a solution 

u =(c 2 M K/2 ,3 /2 (ir 1 ) + c 3 M K/2j _3 /2 (|r 1 ) 

+ (A;/6)(M K/2 ,3 /2 (|r 1 ) 1 M^s/aClr 1 )^ (6.11) 

for 

p= _i jn = -4,Ci = (6.12) 

rea/ constants C 2; C*3, C4. 

7. Group-invariant Radial Solutions 

The radial Schrodinger equation (jl.2p describes a general model for the slow modulation 
of n-dimensional radial waves in a weakly nonlinear, dispersive, isotropic media. Specifically, 
the amplitude a(t,r) of radial harmonic waves with a dispersion relation u(k) is related to 
u(t, r) by p 

a(t, r) tt(ei, e(r — vt)) exp(i(/tr — u;t)), (7.1) 

in which |e| 1 is a small parameter and v = dco/dn is the radial group velocity. 

There is an equivalent planar (i.e. 2-dimensional) form of the radial Schrodinger equation 
(JTSJ) given by 

iu t = u rr + (1 — m)r~ u r + k\u\ p u (7.2) 
containing an extra modulation term mu r /r with a parameter m = 2 — n which is applicable 
for any value of n. The precise interpretation of this extra term can be explained through 
considering the net modulation defined by the 2-dimensional integral of u 

POO 

C(t) = urdr (7.3) 
Jo 

on the whole radial domain in M 2 . If u is regular at the origin (i.e. lim^o u r = 0), then this 
integral formally satisfies 

^ = -(S + \imm), (7.4) 

dt y r^O 

where 

roc 

S(t) = ik u\u\ p rdr (7.5) 
Jo 
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measures the net amount of modulation produced by the nonlinear term in the radial 
Schrodinger equation (17. 2p . Note that C is conserved, — = 0, whenever both the modula- 
tion term and the nonlinear term are absent, m = k = (describing a linear, non- dispersive 
medium). Thus, when m/0, the modulation term mu r /r models the effect of a point-source 
disturbance at the origin r = 0, which alters the slow modulation of radial harmonic waves 
in a planar, weakly nonlinear, dispersive medium. In polar variables u = Aexp(i$), we have 

^ + S = -mlim(Aexp(i($ + 7r/2))), (7.6) 

where C is the net modulation (17.3)) . and also 

dX f°° 

— = -m / 2A 2 $ r dr, (7.7) 
dt J 

where X is the net intensity defined by the integral 

POO 

X(t) = / \u\ 2 rdr. (7.8) 
Jo 

7.1. Solutions from the optimal translation-symmetry subgroup. All of the group- 
invariant solutions derived in section H] using the optimal subgroup of point symmetries 
( I2.41bl) will now be written out in the form u = f(t,r) for the radial Schrodinger equation 
(11. 2p in the next two theorems. Wherever possible, solutions are merged into families that 
do not overlap. 

Theorem 5. The radial Schrodinger equation (11.21) has the following group-invariant solu- 
tions with respect to time-translations ( 12.511) : 

u = (2(p(n - 2) - 2)/kp 2 ) 1/p r ~ 2/p exp(i Cl ), p ^ 2/(n - 2), k/(p(n - 2) - 2) > 0, (7.9) 
u = (±8(p + 2)/(kp 2 )) 1/p (c 2 r 2 ± l/c 2 )- 2 / p exp(i Cl ), 



p = 4/(n - 2), ±k(l - 2/n) > 0, n ^ 2, 

-2/p 



(7.10) 



« =(-|^ 2 (P + 2))-^(r + c 2 r 3 -)-^exp(i Cl ), n 

p = 2(n-3)/(2-n),Jfe/(n-2) <0,n^2,3, '' ' 

" = ( - \{k/n)r 2 + c 3 r 2 ~ n + c 2 ) exp(ici), p = -l,n^0,2, (7.12) 

m = ( — |fcr 2 + c 3 lnr + c 2 ) exp(ici), p=— l,n = 2, (7-13) 

M = (96/A;)(r 2 + c 2 )(2r 2 + c 2 )~ 2 exp(ici), p=l,n = 16, (7.14) 

m = ( — |/cr 2 lnr + c 3 r 2 + c 2 ) exp(ici), p=— l,n = 0, (7-15) 

u =(4A;/3) 1/4 r 2 (r- 2 - c 2 /A;) 3/4 exp (i Cl - i^((fc/c 2 )r- 2 - I)" 1 / 2 ), 
p = -4, n = 0, k > 0, 

it ={-12k) l '\c 2 /k - r 2 / 3 ) 3 / 4 exp (i Cl + i^((£;/c 2 )r 2 / 3 - l)" 1 ^ 
p= -4,n = 4/3, fc < 0, 
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(7.16) 
(7.17) 



((3A;/c 3 2 )r 4/3 - 4c 3 2 (r 1/3 + c 2 r 2 / 3 ) 2 ) 1/4 

x exp ^ici + iiarcsin ((l - |(c 3 4 /A;)(r- 1/3 + c 2 ) 2 ) 1/2 )) , (7-18) 

p = _8,n = 5/3, k > 0, 

(2/A;)r- 2 (2c 2 r 1/3 + 3)(c 2 r 1/3 + 3)- 2 exp(ici), p=l,ra = 13/3, (7.19) 

( - 2c 3 r + (±(* - c, 2 )) 1 / 2 ^ 2 =f l/c 2 )) 1/2 exp(ic 1 ), p = -4,n = l,c 3 2 ^ fc, (7.20) 
( - 2c 4 r + (k - c 3 2 - c 4 2 ) 1 / 2 (c 2 r 2 - l/c 2 )) 1/2 

p = -4, n = 1, c 3 2 + c 4 2 < fc, 

(7.21) 

( - 2c 4 r + (c 3 2 + c 4 2 - A;) 1/2 (c 2 r 2 + l/c 2 )) 1 ' 2 

( in i 2 lN -i/2 i /2(c 3 2 + c 4 2 -/c) 1 / 2 r-c 4 (c 2 r 2 + l/c 2 )\\ 
X 6XP ( 1Cl - ^ - l) ' arCCOSh ( 2c 4 r-(c 32 + c 42 - fe )V 2 ( C2r2 + l/ C2 ) ) J • 



p = -4, n — 1, c 3 2 < k < c 3 2 + c 4 2 , 



(7.22) 

( - 2c 4 r + (c 3 2 + c 4 2 - kfl\c 2 r 2 + l/c 2 )) 1/2 

p = -4 )72 = l,c 3 2 > fc, 

(7.23) 

(-2c 3 r + c 2 r 2 ) 1 / 2 exp(ic 1 ±i|(A;/c 3 2 -l) 1 / 2 ln|c 2 -2c 3 /r|), p = -4,ra = l, (7.24) 
(-2c 3 r + c 2 ) 1/2 exp (i Cl ± i|(A;/c 3 2 - 1) 1/2 In |c 2 - 2c 3 r|), p = -4, n = 1, (7.25) 

(ca/c.) 1 / 2 ^ - l)exp (i Cl - ^(fcV^d^ll)), p = -4,n = 1, * > 0, (7.26) 

c 2 rexp (ici - i(A; 1/2 /c 2 2 )r- 1 ), p = -4, n = I, k > 0, (7.27) 
c 2 exp (ici - i(A; 1/2 /c 2 2 )r), p = -4, n = 1, A; > 0, (7.28) 
(((c 3 2 -fc)/c 4 )r 2 + C4 (l + C2 r) 2 ) 1/2 

x exp ^id - ic 3 (k - c 3 2 )- 1 / 2 arcsinh(( - 1 + (c 4 2 /(A; - c 3 2 ))(r- 1 + c 2 ) 2 )" 1/2 )^, 

p = —4, n = 1, c 3 2 < fc, c 4 > 0, 

(7.29) 
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u=(((c 3 2 -k)/c 4 )r 2 + c,(l + c 2 r) 2 ) 1/2 



u 



x exp [x Cl - ic 3 (c 3 2 - k)- 1 ' 2 arcsin ((l + (c 4 2 /(c 3 2 - k))^ 1 + c 2 ) 2 )" 1/2 )^), 
p = —4, n — 1, c 3 2 > fc, C4 > 0, 

(7.30) 

: (((c3 2 -^)/c4)r 2 + c 4 (l + c 2 r) 2 ) 1/2 



x exp h Cl - ic 3 (k - c 3 2 )- 1/2 arccosh((l + (c 4 2 /(c 3 2 - k))^ 1 + c 2 f) 1/2 ) 
p = —4, n = 1, c 2 < k, c 4 < 0, 

(7.31) 

u = (4(fc - c 3 2 )) V4 (c 2 r 2 ± r) 1 ' 2 exp (i Cl + |ic 3 (/c - c 3 2 )- 1 ' 2 In |c 2 ± r' 1 1) , 

(7.32) 

p= -4,n = l,c 3 < fe, 
m = c 3 (l + c 2 r)exp (i Cl tKA: 172 /^ 2 )^" 1 + C,)- 1 ), p=-4,n = l. (7.33) 

Theorem 6. T/ie radial Schrddinger equation (11.21) /jas t/te following group-invariant solu- 
tions with respect to time-translations combined with phase rotations (12.43bj) : 



u 



u 



(r 1 n/2 (c 2 J| n _ 2 | /2 (v^r) + c 3 ^ n _ 2 | /2 (v^r)) - k/vj exp(i Cl + wt), p = -1, (7.34) 
^r 1_n/2 (c 2 /| n _ 2 | /2 (v /z ^r) + c 3 K\ n _ 2 \/2(V^vr)) - k/vj exp(ici + iut), p = -1. 



(7.35) 

7.2. Solutions from the optimal scaling-symmetry and conformal-symmetry sub- 
groups. All of the group-invariant solutions derived in sections [5] and [6] using the optimal 
subgroups of point symmetries (I2.41cp and (I2.41dj) will be written out in the form u = f(t, r) 



for the radial Schrddinger equation (11. 2p in the next two theorems. 

Theorem 7. The radial Schrddinger equation (II. 2p has the following group-invariant solu- 
tions with respect to scalings combined with phase rotations (I2.43c|) / 



u 



--\ikr 2 (M^ /2 ^ /2 {\w 2 /t) [ ' M^ l%3l2 {\xC x )di 

JC2 

r t/r 2 

- M_ W2j3/2 (iir 2 /t) / M_ lll/2) _ 3/2 (\xC l ) dn exp (i Cl - |i(r 2 /t) + §uilni) , 
p = —1, n = —4. 

(7.36) 
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Theorem 8. The radial Schrddinger equation (jl.2p has the following group-invariant solu- 
tions with respect to time-translations and inversions combined with phase rotations (12.43dj) : 



u 



r-(l+t 2 )/r 2 

#r 2 (M KA3/2 (±r 2 /(l + t 2 )) / M K/2 ^ /2 {\C l )di 



-M K/2 ^ /2 {\r 2 /{l + t 2 )) / M^g/adr 1 )^) (7-37) 

</c 3 

x exp (ici — ik arctan(l/£) — \ir 2 t/ (1 + t 2 )) , 
p = — 1, n = —4. 

7.3. Solutions from the full symmetry group. From Theorem [3l the full group of point 
symmetries admitted by the radial Schrddinger equation (jl.2p can applied to each of the 
solutions u = f(t,r) listed in Theorems [5], El O and Phase rotations (12.31aj) and scal- 
ings (I2.31cp change only the constants appearing in these solutions, while time-translations 
f!2.31bl) at most shift t by a new constant. In contrast, inversions (12.31dp produce additional 
new solutions, which are listed in the next theorem. 

Theorem 9. In the case of the conformal power p = 4/n, the radial Schrddinger equation 
(II. 2p has the following additional group-invariant solutions: 

u =(n(n - 4)/(4A;))"/ 4 r-"/ 2 exp (i Cl - ic 2 r 2 /(4(l + c 2 t))) , 

n ^ 4,Jfe/(n(n-4)) > 0, 

u =(-n 3 /{2k{2n + 4)))"/ 4 r""/ 2 (l + c 2 (r/(l + c 3 t)) 2 - 4/p )~ 2/p exp ( 1Cl - ic 3 r 2 /(4(l + c 3 t))), 
n 2 - n + 4 = 0, 

(7.39) 

u ={\kr 2 + c 3 r 6 (l + c 4 t)" 4 + c 2 (l + c 4 t) 2 ) exp (i Cl - ic 4 r 2 /(l + c 4 t)), 
p = —l,n — —4, 

u =(V 3 (1 + c 4 t)- 1 (c 2 J 3 (v / ^r/(l + CAt)) + c 3 F 3 (v^r/(l + c 4 t))) - (k/u)(l + c 4 tf) 

x exp (ici + u4/(l + c 4 t) - ic 4 r 2 /(4(l + c 4 t))) , ( 7 - 41 ) 
p = — 1, n = —4, i/ > 0, 

« =(r 3 (l + c 4 t)- 1 (c 2 / 3 ( v / ^r/(l + c 4 t)) + c 3 K 3 (^r/(l + c 4 t))) - + c 4 t) 2 ) 

x exp (ici + iz4/(l + c 4 t) - ic 4 r 2 /(4(l + c 4 t))), 
p = — 1, n — —4, v < 0, 



(7.42) 



|ifcr 2 (M_ W2i _ 3/2 (iir 2 /(t(l + c 4 t))) / M_ W2)3 / 2 (K" 1 ) # 

-M_ W2j3/2 (iir 2 /(t(l + c 4 t))) / M.^.^dir 1 )^) (7-43) 

Jc 3 



x exp (ici - |i^ln(c 4 + 1/i) - §ir 2 (l + 2c 4 t)/(t(l + c 4 t))) 
p = —1, n = —4, 
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/ , (t 2 + (1+C4t) 2 )/r 2 

u =\kr 2 {M R / 2 ^ /2 (\r 2 /(t 2 + (1 + c 4 t) 2 )) / M« /2j _ 3 / 2 (|r 1 ) ^ 

(t 2 + (1+C44) 2 )/r 2 

- M K/2i _ 3/2 (|r 2 /(t 2 + (1 + c 4 t) 2 )) / M^^dr 1 ) (7-44) 

x exp (ici — iKarctan(c 4 + 1/t) — ^rr 2 (c 4 + t + c 2 i)/{t 2 + (1 + c 4 t) 2 )) , 
p — —l,n — —4. 

7.4. Analytical behaviour of group-invariant solutions. The radial solutions listed in 
Theorems [5] to [9] exhibit several types of interesting behaviour describing (1) standing waves; 
(2) static and dynamic monopoles; (3) static "bright solitons"; (4) static and dynamic "dark 
solitons" . 

A function u(t,r) describes a radial monopole if it is smooth on < r < oo such that 
\u\ — > as r — > oo and \u\ — > oo as r — > 0. Solution (17. 9p for p > is a n-dimensional static 
monopole; solution (I7.14p for c 2 = is a 16-dimensional static monopole; solutions (17. 9p for 
n ^ 1,2,..., (17. lip for 2 < n < 3, and (17. 19[) are planar static monopoles. Solution (I7.38P 
is a n-dimensional monopole with a dynamic phase; solution (I7.39P is a planar dynamic 
monopole. 

A function u(t, r) describes a radial standing wave if it is smooth on < r < oo such that 
u = U{r) exp(iajt) with u / and \U\ bounded as r — > oo. Solution (I7.34p and solution 
(I7.35P in the case c 2 = are n-dimensional standing waves. 

A function u(t, r) describes a "bright radial soliton" or a "dark radial soliton" if it is smooth 
on < r < oo such that \u\ — > as r — )• oo or \u\ — )• A ^ as r — )■ oo, respectively. Solution 
(17.1 Op for n > 2 in the "+" case is a n-dimensional static bright soliton; solution (I7.14p in 
the case c 2 = is a 16-dimensional static bright soliton. Solution (I7.42p in the case c 2 = 
is a planar dynamic dark soliton exhibiting blow-up. 

More details about these physically interesting radial solutions, including their L 2 norms 
and conserved energies, will be discussed in a separate paper [12] . The remaining radial solu- 
tions found in Theorems E] to O all have unphysical behaviour, in particular \u\ is unbounded 
as r — t- oo. 

8. Concluding Remarks 

In this paper, all explicit group- invariant solutions have been derived (cf Theorems [5] to [9]) 
for the class of semilinear radial Schrodinger equations (II. 2p with a power nonlinearity p ^ 
in multi-dimensions n ^ 1. Among these solutions u(t,r), a few describe n- dimensional 
radial standing waves, radial monopoles, and static radial "bright solitons", which have 
some physical interest. 

Several solutions exist, surprisingly, only for non-integer values of n. In such cases the 
radial Schrodinger equation (II. 2p is shown to have an alternative interpretation as a planar 
(i.e. 2-dimensional) radial equation (II. 2p containing an extra term mu r /r that describes a 
point-source disturbance at the origin r = 0, with m = 2 — n. Some of these planar solutions 
are physically interesting dynamic radial monopoles and dynamic radial "dark solitons" . 

However, no n-dimensional radial solutions are obtained in the analytically important 
case of the conformal power p = 4/n relevant for critical dynamics when n > 2. To look 
for such solutions, we plan to apply the method of group-foliation reduction [13], which 
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has been successfully used in previous work [TJJ [15] to obtain blow-up and dispersive radial 
solutions to semilinear wave equations and semilinear heat conduction equations with power 
nonlinearities in multi-dimensions. 
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